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Abstract 


This  study  evaluates  the  role  of  drift  in  maintaining 
relative  microdistribution  of  larval  populations  of  two 
mayfly  species  in  the  Pembina  River,  Alberta,  Canada. 

Ephemerella  inermis  Eaton  is  univoltine.  Eggs  hatch 
in  August  and  larvae  emerge  the  following  July.  Baetis 
tricaudatus  Dodds  produces  three  cohorts  per  year. 

Emergence  periods  occur  in  early  June,  July  and  early 
September.  Only  one  cohort  overwinters  as  larvae. 

Early  instar  E.  inermis  larvae  are  most  abundant  near 
shore  but  move  to  deeper  water  in  late  autumn.  Baetis 
tricaudatus  larvae  are  seldom  numerous  in  nearshore  areas. 
Microdistribution  of  both  species  is  largely  determined  by 
the  discharge  pattern  of  the  river:  flooding  results  in 
temporary  relocation  to  nearshore  regions. 

Drift  of  both  populations  is  most  marked  in  late  spring 
and  summer.  Nocturnal  drift  densities  are  much  greater  in 
slow  than  in  fast  water  regions. 

Laboratory  studies  indicate  that  departure  from  the 
substrate  is  density-independent  and  independent  of  amount 
of  detritus  in  the  substrate.  Nocturnal  departure  was  more 
pronounced  under  conditions  of  low  mean  current  velocity 
than  in  faster  water.  The  reverse  was  true  during  daylight. 
Rates  of  loss  of  dead  animals  from  the  substrate  were  much 
lower  than  departure  rates  of  live  organisms. 

Settling  rates  of  both  live  and  dead  animals  were 
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inversely  related  to  mean  current  velocity.  Live  animals 
had  greater  settling  capacity  than  did  dead  organisms. 

Lateral  dispersal  of  suspended  larvae  depended  on  both  the 
downstream  distance  travelled  and  on  the  distance  from  shore 
that  larvae  entered  the  water  column.  Live  animals  became 
more  markedly  dispersed  laterally  than  did  dead  animals. 

A  model  was  developed  to  predict  changes  in  benthic 
distribution  of  the  two  populations  that  result  from  drift. 
Lateral  transport  of  drifting  animals  was  assumed  to  be  a 
Markov  process.  Probability  estimates  of  departure,  down¬ 
stream  movement  and  lateral  dispersal  were  combined  to 
produce  transition  matrices.  Temporal  changes  in  micro¬ 
distribution  across  the  river  were  then  estimated  as  a 
Markov  chain.  Although  predicted  estimates  of  drift  broadly 
corresponded  with  observed  levels,  predicted  and  observed 
benthic  distributions  were  markedly  different.  I  conclude 
that  drift  should  result  in  rapid  transport  of  most  individuals 
in  each  population  to  stream  margins. 

I  suggest  that  drift  is  an  important  determinant  of 
microdistribution,  but  observed  distributions  of  larvae 
are  maintained  only  by  the  ability  of  animals  to  crawl 
towards  the  river  centre.  Active  drift  of  B.  tricaudatus 
and  E.  inermis  larvae  may  be  an  adjunct  to  normal  crawling 
search  behaviour.  Drift  exhibited  by  relatively  sedentary 
taxa  must  serve  other  functions. 
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INTRODUCTION 


Intensive  research  on  aspects  of  invertebrate 
drift  in  recent  years  has  done  much  to  elucidate  factors 
influencing  downstream  transport  of  benthic  animals. 
Representatives  of  all  species  can  occasionally  be  found 
suspended  in  the  water  column.  This  may  be  the  incidental 
result  of  catastrophic  changes  in  the  river  (catastrophic 
drift).  Alternatively,  drift  may  appear  in  recurring 
seasonal  or  diel  cycles,  independently  of  visible  mod¬ 
ifications  in  streamflow  characteristics.  Occurrence  of 
the  latter  is  referred  to  as  behavioural  drift  and  is  most 
prevalent  in  herbivorous,  hemimetabolous  aquatic  insects, 
notably  Ephemeroptera.  The  prevailing  view  of  drift 
attributes  these  periodicities  to  an  interaction  between 
cyclic  nocturnal  or  seasonal  increases  in  the  activity  of 
individuals  and  the  eroding  influence  of  the  water  current. 
This  causes  an  accidental  dislodgement  of  animals. 

Interpretations  of  the  adaptive  significance  of  the 
phenomenon  differ.  Researchers  who  have  studied  streams 
in  which  the  incidence  of  drift  is  low  relative  to  benthic 
densities  ascribe  little  ecological  importance  to  drift. 
However,  in  many  rivers  the  number  of  animals  drifting 
over  a  portion  of  substrate  during  24  hours  frequently 
exceeds  the  benthic  density.  Muller  (1954)  viewed  drift 
as  a  means  of  reducing  competition  for  food  or  space  when 
high  densities  of  animals  occur  on  the  stream  bottom. 
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Individuals  that  enter  the  drift  then  become  available  to 
colonize  unexploited  areas  downstream.  Drift  may  thus 
serve  the  dual  role  of  regulating  and  dispersing  pop¬ 
ulations.  Muller  (195^)  felt  that  some  compensating 
mechanism  was  required  to  account  for  an  attrition  of 
populations  in  upstream  areas  and  proposed  a  "colonization 
cycle",  wherein  emerging  adults  flew  upstream  to  oviposit 
and  repopulated  these  reaches. 

Waters  (I96I)  observed  that  drift  in  highly  productive 
rivers  was  much  greater  than  in  sparsely  populated  ones, 
but  upstream  areas  seldom  became  depopulated  by  drift. 
Consequently,  he  felt  upstream  flight  was  not  a  requisite. 
Waters  postulated  the  alternative  hypothesis  that  drift 
in  response  to  competition  serves  primarily  to  remove  those 
individuals  which  the  habitat  cannot  support.  Drifting 
animals  are  thus  indicative  of  excess  production. 

Another  view  holds  that  drift  represents  the  net 
response  of  individuals  to  the  suitability  of  their 
respective  microhabitats.  Early  field  studies  revealed 
an  exponential  relationship  between  total  invertebrate  drift 
and  stream  discharge  or  the  strength  of  the  current  (Bailey 
1966,  Elliott  1967b).  However,  other  workers  have  found 
that  reduced  current  can  also  produce  increased  drift 
(Edington  1965>  Elliott  1968a,  Minshall  and  Winger  1968, 
Corkum  et  al .  1977,  Gore  1977).  Such  evidence  suggests 
that  departure  from  the  substrate  may  be  active  at  times 
rather  than  accidental,  but  more  signif icantly  it  implies 
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that  drift  is  a  controlled  response  to  an  environmental 
gradient.  The  range  of  current  velocity  within  which 
minimal  drift  occurs  may  be  representative  of  the  zone 
of  preference  of  a  particular  population;  outside  this 
zone  animals  tend  purposefully  to  leave  their  microhabitat. 

Similar  observations  have  been  made  on  the  influence 
of  other  abiotic  factors  affecting  drift,  e.g.  substrate 
texture,  water  temperature  and  oxygen  concentration.  The 
influence  of  predators  (Corkum  and  Pointing  1979 »  Corkum 
and  Clifford  1980,  Peckarsky  1980b,  Walton  1980b)  and 
food  availability  (Hildebrand  197^ »  Keller  1978,  Bohle 
1978,  Walton  1978)  have  also  been  implicated  in  modifying 
the  magnitude  of  drift.  Thus,  under  certain  circumstances, 
benthic  organisms  may  drift  in  response  to  many  of  the  biotic 
and  abiotic  variables  typically  used  to  define  aquatic 
microhabitats . 

Figure  1  summarizes  the  factors  that  have  been  shown 
to  influence  drift.  Pathways  are  numbered  and  corresponding 
references  to  them  are  given  in  Table  1.  The  tabular  summary 
lists  only  those  studies  in  which  a  'positive'  relationship 
was  found  or  inferred.  General  descriptive  accounts  of 
field  studies  have  not  been  included. 

The  relative  importance  of  each  of  the  parameters 
illustrated  in  Fig.  1  has  yet  to  be  determined  for  any 
species,  but  it  is  unlikely  that  all  are  prominent  in 
governing  the  drift  of  a  particular  population.  The  complexi 
of  the  figure  underlines  the  central  problem  inherent  in  a 


X  ■  1,1  £Bmirf-to 

,  ' 


- 


' 

■ 


4 


d 

s 

3 

I — 1 


0 

d 

0 

& 


o  e  s  cd 
O  o  3  4° 

H  Ch  O 
0 

•H  T)  -Pi 
Cd  cd  CD  O 

^00  »rH 


0 

Pi  CD 

+"  £ 
o 
Sh  d 
0  u 

d  ^ 

0 


0  4ft 

d  £ 

o 

P-t  tH 
0 

0  0 
d  |H 
P 
d  cd 
cd  Eri 


0  0  3d 
H  P  O  .H 

cd  cd 

S  d  p  0 

•H  P  d  0 

d  0  cd  o 
cdp  a  d 
3  0 

4ft  0  0  d 
O  o  0 
•H  0  3  <P 
0 

£  -P  d  d 
•H 

>5  0  o 
4ft  P  P  -P 
cd 

0  d  P  ^ 

£  d  p  £ 

0  3  O 
•H  P  «H  P< 
d  0  H  0 
Cd  d  3  0 
4ft  £  d 
Od.H  d 
0  d  p  o 
E  cd  0  o 


0 

d 

3 

W) 

•H 


the  indicated  pathways. 


Table  1.  Summary  of  studies  on  factors  that  influence 

departure  from  substrate  or  downstream  transport 
of  aquatic  invertebrates.  Numbers  heading  each 
column  refer  to  pathways  shown  in  Fig.  1. 
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complete  understanding  of  the  initiation  of  drift. 

In  general,  the  importance  of  active  entry  into  the 
water  column  has  probably  been  understated  in  much  of  the 
drift  literature  (Walton  1980a) ,  but  whether  entry  occurs 
by  either  active  or  by  passive  means,  mediated  by  activity 
changes,  is  of  limited  importance  to  the  outcome  of  the  act. 
Ulf strand  (I968)  stated  that,  regardless  of  the  means  by 
which  animals  become  suspended  in  the  water  column,  ’’the 
start  of  drifting  movement  is  by  no  means  accidental." 
Townsend  and  Hildrew  (1976)  observed  that  the  benefits  of 
drift  to  an  individual  will  accrue  regardless  of  whether 
entry  into  the  water  column  is  accidental  or  intentional. 

Once  in  the  water  column,  organisms  are  carried  down¬ 
stream.  A  constant  proportion  return  to  the  bottom  for 
every  unit  of  distance  or  time  travelled  (McLay  1970 , 

Elliott  1971a,  McKone  1975)  •  Reestablishment  may  be 
facilitated  by  swimming  motions  made  by  the  animals  (McKone 
1975).  Alternatively,  downstream  transport  of  certain 
animals  may  be  passive  (Ulf strand  1968)  and  terminated  by 
thigmotactic  responses  of  the  individuals  upon  chance  contact 
with  the  substrate.  A  behavioural  component  of  settling 
results  in  increased  rates  of  return  to  the  bottom  by 
certain  species  (Elliott  1971a,  Walton  1978,  Ciborowski  and 
Corkum  I98O) .  Animals  able  to  stage  a  rapid  return  to  the 
substrate  tend  to  be  those  that  frequently  show  diel  drift 
periodicities  (Ciborowski  and  Corkum  I98O) . 

Current  exerts  the  most  important  influence  on  settling 
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rates,  although  substrate  texture,  water  temperature,  benthic 
density  at  site  of  reestablishment,  and  the  presence  or 
absence  of  light  may  be  modifying  factors  (references  in 
Table  1 ) . 

A  fraction  of  drifting  animals  are  lost  as  a  result  of 
predation  by  fishes,  emergence,  and  mechanical  injury. 

Mancini  et  al .  (1979)  found  that  removal  of  drift-feeding 
fishes  from  a  pool  had  no  influence  on  the  density  of 
invertebrates  drifting  through  it.  Consequently,  although 
drifting  animals  may  be  an  important  food  source  for  fish 
predators,  the  probability  of  an  individual  insect  being 
eaten  is  low.  Laboratory  studies  suggest  that  mortality 
due  to  mechanical  injury  is  minimal,  even  at  very  high 
water  velocities  (personal  observations).  Thus,  most 
individuals  that  leave  the  substrate  probably  return 
successfully. 

Drift  often  provides  the  primary  source  of  animals 
during  recolonization  of  disturbed  or  denuded  stream  areas 
(Waters  1964,  Townsend  and  Hildrew  1976,  Williams  and  Hynes 
1976,  Sheldon  1977,  Gore  1979,  Wise  and  Molles  1979,  Shaw 
and  Minshall  I98O) .  When  adjacent  upstream  reaches  support 
high  densities  of  individuals,  normal  standing  stocks  of 
animals  can  be  accumulated  on  a  scoured  substrate  in  as 
little  as  4  days  (Waters  1964) .  Such  recolonization  must 
also  be  a  reflection  of  the  turnover  of  individuals  that 
occurs  in  all  parts  of  a  river  by  drift  and  subsequent 
resettling. 

The  number  of  animals  on  a  portion  of  river  bottom 
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necessarily  represents  the  balance  between  net  immigration 
and  emigration  as  well  as  births  and  deaths.  Since  relatively 
minor  changes  in  stream  discharge  can  produce  large  changes 
in  the  flow  pattern  over  the  bottom  and  consequently  influence 
both  the  substrate  composition  and  food  resources,  it  may 
be  expected  that  animals  are  adapted  to  relocate  when 
conditions  become  unfavourable.  Drift  could  be  an  important 
means  by  which  this  relocation  occurs. 

My  study  was  undertaken  to  evaluate  the  importance  of 
drift  in  governing  the  microdistribution  of  populations  of 
two  mayfly  species,  Baetis  tricaudatus  Dodds  and  Ephemerella 
inermis  Eaton,  in  the  Pembina  River,  central  Alberta,  Canada. 
Both  species  are  numerically  abundant  in  streams  throughout 
western  North  America  and  frequently  dominate  the  drift  in 
these  streams. 

I  have  presented  my  research  in  four  chapters.  The 
first  establishes  the  life  history,  microdistribution  and 
drift  characteristics  of  larvae  at  a  site  on  the  Pembina 
River  near  Entwistle,  Alberta.  Chapter  2  reports  the 
results  of  laboratory  studies  which  assess  the  relative 
importance  of  current  velocity,  benthic  density  of  larvae 
and  detritus  concentration  on  departure  rates  of  the  organisms. 
I  also  studied  the  effects  of  current  velocity  on  washing  away 
of  dead  animals  to  examine  the  role  that  behaviour  might 
play  in  inducing  animals  to  leave  the  substrate.  Chapter 
3  provides  information  on  the  transport  of  larvae  in  the 
Pembina  River.  Live  and  dead  animals. were  released  into 
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the  water  column.  Catches  of  larvae  in  a  series  of  down¬ 
stream  nets  permitted  me  to  estimate  rates  of  settling  and 
the  degree  of  lateral  transport  that  drifting  animals  are 
subject  to.  In  Chapter  4,  I  derive  a  mathematical  model 
which  incorporates  data  on  larval  microdistribution, 
departure  rates,  downstream  transport,  and  lateral  transport 
to  examine  how  drift  may  modify  the  microdistribution  of  the 
two  larval  populations.  Predictions  of  the  relative  distrib¬ 
ution  of  larvae  after  2  weeks'  drift  and  the  number  of  animals 
expected  in  the  drift  were  compared  with  data  collected 
from  the  Pembina  River. 
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CHAPTER  1 


LIFE  HISTORIES,  MICRODISTRIBUTION  AND  DRIFT 
OF  TWO  MAYFLY  (EPHEMEROPTERA)  SPECIES 
IN  THE  PEMBINA  RIVER, 

ALBERTA,  CANADA. 
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Introduction 

Baetis  tricaudatus  Dodds  and  Ephemerella  inermis 
Eaton  are  among  the  most  widespread  and  abundant  mayfly 
species  in  western  North  America  (Johnson  1978,  Morihara 
and  McCafferty  1979)*  Although  E.  inermis  has  a  simple 
univoltine  life  cycle  (Hartland-Rowe  1964,  Zelt  1970* 

Hamilton  1979*  Barton  1980),  its  microdistribution  may 
vary  with  its  developmental  stage  (Minshall  and  Minshall 
1977*  Rabeni  and  Minshall  1977) •  Corkum  and  Pointing 
(1979)  reviewed  development  of  B.  tricaudatus  (as  B.  vagans 
McD.)  in  eastern  North  America,  but  the  life  cycle  of  western 
populations  is  complex  and  poorly  understood  (Clifford  et  al. 
1973) .  Larvae  of  this  species  are  normally  found  in  fast 
water  (Minshall  and  Minshall  1977,  Rabeni  and  Minshall  1977* 
Corkum  and  Pointing  1979) • 

Both  E.  inermis  and  B.  tricaudatus  larvae  are  frequently 
cited  as  dominant  organisms  in  stream  communities,  yet 
studies  on  their  biology  usually  treat  only  their  life  histories 
or  their  microdistribution,  seldom  both.  I  investigated 
whether  microdistribution  of  these  two  mayflies  changes  as 
the  larvae  develop;  and,  if  so,  what  mechanisms  are  responsible 
for  relocation  of  the  animals. 


7  *;  4^ 


' 


' 


Methods 


Study  Site 

The  study  was  conducted  on  the  Pembina  River  near  the 
town  of  Entwistle,  central  Alberta,  Canada  (53°37  N., 

115°00  W).  This  fifth  order  river  originates  in  the  Rocky 
mountains  and  flows  into  the  Athabasca  River.  At  the  study 

p 

site  it  drains  a  watershed  area  of  4530  km  ,  is  65  m  wide, 
and  has  a  maximum  depth  of  1.2  m  during  normal  flows. 
Substrate  is  uniform  across  the  channel,  consisting  of  smooth 
limestone  cobble  embedded  in  coarse  sand.  The  substrate  is 
barren  of  macrophytes  except  during  extended  periods  of  low 
discharge,  at  which  time  Cladophora  sp.  may  become  abundant. 
Current  velocity  increases  evenly  from  shore  to  stream 
centre.  Cross-sectional  mean  velocity  is  0.4  m  s  during 
typical  summer  water  levels.  Bankside  vegetation  (Populus 
balsamif era  Li ,  Alnus  rugosa  (Du  Roi)  ,  Salix  spp .  ,  Picea 
spp. )  provides  the  majority  of  detrital  material.  The  river 
is  ice-covered  from  November  until  April.  Kellerhals  et  al . 
(1972)  summarized  the  hydraulic  and  geomorphic  characteristic 
of  the  Pembina  River  above  the  study  site. 

Benthos  Collection 

Animals  were  collected  with  a  modified  Hess  sampler 

2 

(Waters  and  Knapp  I96I)  that  enclosed  a  sample  area  of  0.1  m 
Pore  size  of  the  netting  was  0.380  mm.  Benthic  samples  were 
taken  at  10-15  day  intervals  during  the  ice-free  period  from 
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May  1978  to  March  1980.  Collections  were  made  along  each 
of  three  transects  running  parallel  to  the  shoreline. 

Transects  on  each  sampling  date  were  chosen  on  the  basis 
of  water  depth;  at  10-15  cm,  25-35  cm  and  45-55  cm.  Except 
during  periods  of  extremely  high  or  low  flow,  these  corresponded 
to  distances  from  shore  of  3»  7>  and  11  m.  Current  velocity 
was  measured  at  each  transect  with  a  Gurley,  Stevens  or 
Ott  C-l  current  meter  except  when  heavy  flood  conditions 
released  enough  detritus  to  clog  the  propeller.  Three 
winter  collections  were  taken  through  the  ice,  25  and  32  m 
from  shore.  Water  depth  beneath  the  ice  was  45  and  55  cm, 
respectively,  in  December,  but  had  decreased  to  5  and  10  cm 
by  March. 

Three  Hess  samples  were  taken  1  m  apart  and  pooled  at 
each  transect  in  1978.  In  1979 »  five  separate  samples  were 
collected  along  each  transect.  The  intermediate  transect 
was  not  sampled  in  May  1979*  All  collections  were  preserved 
in  Kahle's  fluid  and  sorted  in  the  laboratory. 

In  1979  1  all  non-animal  organic  material  was  removed 
from  the  inorganic  fraction  by  elutriation,  separated  into 
coarse  (>1.00  mm)  and  fine  (<1.00  mm)  portions,  air-dried, 
and  weighed  to  the  nearest  0.01  g. 

Water  temperature  was  measured  with  a  hand-held  ther¬ 
mometer  between  1100  and  1130  h  on  each  sampling  occasion. 

Daily  discharge  data  were  available  from  a  gauging  station 
maintained  by  the  Water  Survey  of  Canada,  1  km  upstream  of 
the  sampling  site. 
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Life  History 

Size  distribution  of  animals  on  each  sampling  date  was 
determined  by  measuring  head  width  at  the  level  of  the  eyes. 

In  1978 »  up  to  50  animals  per  species  were  measured  from 
each  of  nearshore  and  remote  transects.  Up  to  20  animals 
per  sample  from  all  transects  were  measured  in  1979.  Where 
subsampling  of  collections  was  necessary  for  measurements, 
individuals  were  visually  ranked  by  size  and  the  appropriate 
number  of  organisms  selected  at  regular  size  increments. 

Head  widths  were  measured  with  an  ocular  micrometer  under 
a  dissecting  microscope  at  40X  magnification  to  the  nearest 
0.024  mm. 

Periods  of  emergence  were  estimated  by  sweeping  bank- 
side  vegetation  for  adults  and  by  noting  the  presence  of 
mature  larvae  (those  with  black  wing  pads)  in  benthic  samples. 

Microdistribution 

One-way  analysis  of  variance  (ANOYA)  was  used  to  test 
for  differences  in  density  of  each  species  between  transects 
on  each  date  in  1979*  Where  there  were  significant  differences, 
effects  were  partitioned  using  Student-Newman-Keuls  tests 
(Sokal  and  Rohlf  1969). 

Least-squares  regression  analysis  was  used  to  determine 
the  relationship  between  amounts  of  organic  material  and 
benthic  density  of  each  species.  Separate  regressions  were 
conducted  for  large  (CPOM)  and  small  (FPOM)  particle  sizes. 

In  addition,  the  relationship  between  numbers  of  all  other 
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benthos  (AOB)  in  each  sample  and  the  two  mayfly  species  was 
tested.  The  regression  analyses  were  performed  on  the  data 
pooled  from  all  1979  samples,  on  samples  pooled  by  date,  on 
samples  pooled  by  transect  and  on  sets  of  replicates  collected 
along  each  transect  on  each  date.  Pooled  coefficients  of 
determination  were  calculated  for  each  grouping  of  data  by 
adding  together  all  explained  sums  of  squares  for  a  series 
of  regressions  and  dividing  by  the  pooled  total  sums  of 
squares.  This  provided  an  estimate  of  the  proportion  of 
the  overall  variation  removed  by  regressions  on  the  various 
data  groupings. 

Drift 

Twenty-four  hour  drift  studies  were  carried  out  June 
5-6,  July  5-6  and  August  2-3.  1978.  Drift  was  monitored 
continuously  by  means  of  nets  with  a  mouth  opening  of  9  x  9 
cm,  and  net  pore  size  of  0.320  mm  (Clifford  1972a).  On 
each  occasion  two  nets  were  anchored  in  the  river  with 
steel  rods  driven  into  the  riverbed;  one  4  m  from  shore  in 
15-20  cm  of  water  and  the  other  12  m  from  shore  at  a  depth 
of  45-55  cm.  Each  net  was  kept  3  cm  below  the  water  surface 
and  was  emptied  at  30-min  intervals.  Water  velocity  at  the 
mouth  of  each  net  was  measured  every  6  h.  In  July  and  August 
a  third  net  was  anchored  3  cm  from  the  substrate  on  the  same 
stakes  as  the  deepwater  net. 

Net  contents  were  preserved  in  Kahle's  fluid  and  sorted 
in  the  laboratory.  Because  the  amount  of  water  passing 
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through  the  nets  varied  with  net  position  and  date,  drift 

was  expressed  as  drift  density  (no.  animals  m  ^  water  filtered). 

Results 

Physical  Characteristics 

Daily  discharge  and  water  temperature  differed  sub¬ 
stantially  between  1978  and  1979  (Fig.  2).  Mean  annual  daily 
discharge  in  1978  was  the  second  highest  since  1958;  whereas 
in  1979  it  was  the  second  lowest  (Water  Survey  of  Canada, 
unpubl.  data).  In  1979 >  there  was  relatively  little  var¬ 
iation  in  flow  apart  from  the  normal  run-off  peaks  of  May 
and  June.  Four  additional  major  flood  peaks  occurred  in  1978. 

The  Pembina  River  remained  ice-covered  2  weeks  longer 
in  1979  than  in  1978.  This  delayed  the  spring  warming  of 
the  river.  However,  the  frequent  flooding  in  1978  resulted 
in  lower  summer  water  temperatures  than  in  1979* 

Table  2  summarizes  mean  current  velocity,  detrital 
concentrations  and  benthic  densities  of  animals  other  than 
B  tricaudatus  and  E.  inermis  along  the  three  transects 
sampled  in  1979.  Current  velocities  varied  little  through¬ 
out  the  year  (a  consequence  of  the  sampling  design)  but  FPOM, 
CPOM  and  AOB  were  substantially  greater  during  late  summer 
and  early  autumn  than  at  other  times. 

Life  History 

Ephemerella  inermis  was  univoltine  in  the  Pembina  River. 


. 

* 


Figure  2.  Mean  daily  discharge  (top)  and  water  temperature 


(bottom)  at  sample  site  on  Pembina  River,  1978 
(solid  line)  and  1979  (broken  line). 
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Table  2.  Means  ±  1  standard  deviation  of  current  velocities,  detrital  fractions 
and  density  of  all  other  benthos  along  three  transects. 
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Early  instar  larvae  first  appeared  at  the  end  of  July  1979, 
but  not  until  mid-August  1978  (Fig.  3) •  Maximum  densities 
of  animals  were  reached  in  late  October  1979  and  in  early 
October  1978  (Fig.  4)  .  Growth  had  ceased  by  the  end  of 
October  and  did  not  resume  until  ice  break-up.  Mature 
larvae  first  appeared  at  the  beginning  of  June  1978  but  were 
delayed  10  days  in  1979-  Peak  emergence  occurred  at  the  end 
of  June,  although  some  larvae  persisted  in  the  river  until 
the  end  of  August. 

Larvae  of  B.  tricaudatus  were  present  throughout  the 
year,  although  very  few  remained  in  the  river  in  June  1979 
(Fig.  5) •  Low  water  temperatures  in  early  1979  appeared  to 
delay  hatching  of  eggs  in  spring  more  than  it  delayed  emer¬ 
gence  of  adults.  Major  periods  of  emergence  were  observed 
at  the  end  of  May  and  August.  Lower  numbers  of  individuals 
also  emerged  throughout  July  (Fig.  5) •  The  hatching  of 
young  larvae  was  marked  by  subsequent  peaks  of  benthic  density 
in  June,  August  and  October  (Fig.  6). 

Microdistribution 

During  1978,  spring  and  summer  densities  of  E.  inermis 
larvae  in  slow  water  near  shore  were  as  great  as  or  greater 
than  at  either  of  the  more  remote  transects  (Fig.  4)  .  Near¬ 
shore  densities,  however,  fell  noticeably  during  the  one 
interflood  sampling  period.  Regions  of  greatest  density 
were  reversed  (i.e.  greatest  along  the  remote  transect) 
during  the  same  period  in  1979-  Early  instar  larvae  were 
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Figure  3* 


Frequency  distribution  of  head  widths  of  E.  inermis 


larvae  in  1978  (upper)  and  1979  (lower), 
line  above  histogram  indicates  abundant 
broken  line  indicates  sporadic  emergence 
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emergence , 
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Figure  4.  Benthic  density  of  E.  inermis  larvae  in  1978  (upper) 
and  1979  (lower;  mean  ±  1  standard  error).  Solid 
line  =  nearshore,  broken  line  =  intermediate,  dotted1 
line  =  remote  transect.  Vertical  bars  at  top  of 

j 

lower  figure  connect  transects  for  which  mean 

density  did  not  significantly  differ  (S-N-K 

test,  p  >  0.05)  on  each  1979  sampling  date; 

n,  i,  and  r  represent  nearshore,  intermediate 

and  remote  transects,  respectively.  j 
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Figure  5. 


Frequency  distribution  of  head  widths  of  B. 


tricaudatus  larvae  in  1978  (upper)  and  1979  (lower). 
Remaining  explanation  as  in  Fig.  3* 
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Figure  6.  Benthic  density  of  B.  tricaudatus 
(upper)  and  1979  (lower;  mean  ±  1 
Remaining  explanation  as  in  Fig. 


larvae  in  1978 
standard  error) . 
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first  found  at  slow-water  and  intermediate  transects  and  only 
later  at  the  fast-water  transect.  Appearance  along  the  remote 
transect  was  accompanied  by  a  drop  in  densities  near  shore. 

Spring  and  summer  distributional  patterns  of  B.  tricaudatus 
larvae  were  similar  to  those  of  E.  inermis  (Fig.  6).  Animals 
were  equally  abundant  at  all  transects  in  1978,  but  significantly 
greater  densities  were  found  in  fast  water  along  the  remote 
transect  during  1979*  Densities  at  the  nearshore  transect  did 
not  exceed  those  in  faster  water  as  was  observed  for  E.  inermis 
in  autumn. 

Regression  analysis  of  the  relationship  between  the  three 
independent  variables  (FPOM,  CPOM,  AOB)  and  E.  inermis/B . 
tricaudatus  abundance  removed  the  greatest  amount  of  variation 
when  transect  replicates  within  each  date  were  separately 
analysed  and  then  pooled  (Table  3) •  Positive  relationships 
were  found  between  all  three  variables  and  density  of  E.  inermis 
larvae,  but  only  CPOM  and  AOB  removed  a  significant  amount  of 
the  variation  in  B.  tricaudatus  densities.  The  relationship 
between  the  two  detrital  fractions  and  E.  inermis  abundance 
was  not  evident  when  other  groupings  of  the  data  were  analysed. 

In  contrast,  relationships  between  all  three  independent  var¬ 
iables  and  B.  tricaudatus  larvae  were  negative.  The  partitioning, 
by  transect,  of  the  two  mayfly  species  into  regions  of  similar 
current  velocity  (earlier  1-way  ANOYA  tests)  accounted  for  42 
and  67  percent  of  between  replicate  variation  for  E.  inermis 
and  B.  tricaudatus  larvae,  respectively. 
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Table  3.  Pooled  coefficients  of  determination  for  regressions  of  three  variables 
on  densities  of  B.  tricaudatus  and  E.  inermis  larvae.  *  indicates  that 
regressions  remove  a  significant  (p<.OlT  amount  of  variation  by  ANOVA. 
Symbols  in  parentheses  indicate  whether  significant  regression  coefficients 
are  positive  or  negative. 
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Drift 

Both  E.  inermis  and  B.  tricaudatus  larvae  showed  nocturnal 
increases  in  drift,  especially  in  nearshore  samples.  Although 
daytime  drift  density  of  each  species  was  comparable  at  both 
nearshore  and  remote  sites  for  any  particular  date,  nocturnal 
drift  densities  were  substantially  greater  near  shore  than  in 
deeper  water  (Fig.  7) •  Drift  densities  of  E.  inermis  larvae 
were  greater  than  those  of  B.  tricaudatus  in  June,  the  only 
month  when  benthic  densities  of  the  two  species  were  comparable. 
Drift  densities  of  both  species  were  exceptionally  high  when 
compared  with  values  in  other  field  studies  (see  Armitage  1977)* 

Ephemerella.  inermis  larvae  displayed  a  single  nocturnal 
peak  in  drift  density  in  June,  whereas  B.  tricaudatus  exhibited 
at  least  two  peaks,  the  major  one  occurring  just  prior  to  dawn 
(Fig.  8) .  Patterns  of  drift  on  other  dates  were  similar. 

Discussion 

In  the  Pembina  River,  eggs  of  E.  inermis  hatch  between 
the  end  of  July  and  mid-August.  Larvae  grow  rapidly  in  autumn, 
complete  development  the  following  spring,  and  emerge  at  the 
beginning  of  July.  Similar  developmental  patterns  have  been 
found  in  Alberta  rivers  by  Zelt  (1970),  Hamilton  (1979)  and 
Barton  (1980).  Hartland-Rowe  (1964)  studied  the  life  history 
of  E.  inermis  in  a  spring-fed  foothills  stream  and  found 
emergence  occurring  several  weeks  later. 

Baetis  tricaudatus  was  characterised  by  three  peaks  of 


.  • 

' 

' 

1 


Figure  7*  Drift  density  of  E.  inermis  (upper)  and  B.  tricaudatus 

(lower)  larvae  in  slow  and  fast  water.  U  and  L 
indicate  samples  of  upper  and  lower  portions  of 
water  column.  Histograms  within  each  group  rep¬ 
resent  drift  density  calculated  over  8-h  intervals; 
left,  1300-2100  h;  centre,  2100-0500  h;  right, 
0500-1300  h. 
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Figure  8. 


Drift  density  of  E.  inermis  (broken  line)  and  B. 
tricaudatus  (solid  line)  larvae  in  nearshore  net, 
June  5-6  1978.  Horizontal  bar  represents  darkness. 
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abundance  (November,  August,  late  June)  followed  by  two  major 
(May,  late  August/early  September)  and  a  minor  (July)  period 
of  emergence.  This  species  may  be  trivoltine  in  the  Pembina 
River.  However,  the  appearance  of  young  larvae  concurrently 
with  adult  emergence  in  spring  1978  and  autumn  of  both  years 
brings  this  interpretation  into  question.  Murphy  (1922)  found 
that  B.  tricaudatus  eggs  in  a  New  York  stream  required  28  days 
to  hatch  during  May.  In  the  laboratory,  eclosion  occurred  in 
1 1  days . 

It  is  equally  possible  that  the  life  history  pattern  of 
B.  tricaudatus  is  similar-  to  that  found  by  Murphy  (1922). 

Murphy  reported  that  the  New  York  population  completed  three 
generations  in  two  years  (Fig.  9a).  Two  of  the  cohorts,  which 
hatched  in  autumn,  required  9  months  to  develop  from  egg  to 
adult,  whereas  eggs  laid  in  the  spring  required  only  6  months. 
Clifford  et  al.  (1973)  have  observed  that  mayfly  species  tend 
more  towards  summer  life  cycles  in  subarctic  regions  than  in 
more  temperate  areas.  They  reported  winter  dormancy  of  B. 
tricaudatus  eggs  in  3  Alberta  streams.  It  is  thus  possible 
that  eggs  laid  in  autumn,  which  would  normally  hatch  in 
November  in  more  southern  localities,  overwinter  as  eggs  in 
the  Pembina  River  (Fig.  9b).  Larvae  that  appear  in  June  grow 
rapidly  to  emerge  through  July.  The  offspring  of  this  cohort 
might  not  hatch  until  autumn,  and  then  overwinter  as  larvae. 

The  long  winter  and  rapid  warming  of  water  in  spring  would  tend 
to  produce  a  relatively  synchronized  emergence  of  adults  once 
conditions  became  favourable.  Eggs  of  this  winter  cohort  would 
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Figure  9.  A.  Development  of  B.  tricaudatus  in  Ithaca,  N.Y. 

(after  Murphy  1922). 

B.  Possible  developmental  pattern  of  B.  tricaudatus 
in  Pembina  River.  Stippled  areas  represent  eggs; 
solid  sections  represent  larvae;  arrows  denote 
adult  emergence  and  oviposition. 


38 


39 


then  hatch  in  July  to  produce  another  fairly  strongly  synchronized 
emergence  period  in  late  August/early  September. 

Neither  the  trivoltine  nor  the  alternative  interpretation 
of  B.  tricaudatus  life  history  in  the  Pembina  River  is  consistent 
with  the  egg  development  times  given  by  Murphy  (1922).  If  B. 
tricaudatus  is  trivoltine,  then  I  would  estimate  that  eggs 
would  have  to  hatch  in  less  than  5-7  days  during  the  summer, 
whereas  the  alternative  interpretation  would  require  egg  de¬ 
velopment  times  of  approximately  35  days.  Controlled  studies 
of  egg  and  larval  development  will  be  required  before  the  life 
history  of  this  species  can  be  properly  evaluated. 

Several  workers  have  described  the  life  history  of  B. 
tricaudatus  as  bivoltine  (Robertson  1967 »  Clifford  1969  * 

Zelt  1970,  Hamilton  1979) »  but  all  had  examined  material  collected 
at  3  to  4-week  intervals  from  qualitative  benthic  samples. 

Although  their  data  were  not  precise  enough  to  permit  them  to 
distinguish  the  rapidly  developing  summer  cohort  larvae  from 
the  other  two,  their  data  are  consistent  with  those  from  the 
Pembina  River. 

Microdistribution  of  both  species  reflected  the  different 
flow  conditions  during  the  2  years  of  the  study.  In  1978,  a 
year  of  frequent  heavy  flooding  during  early  summer ,  animals 
were  equally  abundant  at  all  transects.  Larvae  were  much  less 
numerous  near  shore  the  following  year,  when  flow  conditions 
were  more  stable.  B.  tricaudatus  larvae  are  most  commonly 
found  in  fast  water  areas  (Minshall  and  Minshall  1977 »  Rabeni 
and  Minshall  1977,  Corkum  and  Pointing  1979,  Morihara  and 
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McCafferty  1979)*  However,  Lehmkuhl  and  Anderson  (1972)  observed 
that  heavy  winter  flooding  in  an  Oregon  stream  caused  the 
drifting  of  B.  tr i cauda tus  larvae  from  riffle  areas  to  back¬ 
waters.  The  return  of  stable  flow  in  the  spring  resulted  in 
a  movement  of  animals  back  to  riffles. 

Microhabitat  preferences  of  E.  inermi s  were  not  as  distinct 
as  for  B.  tricaudatus  larvae.  In  1979,  most  E.  inermis 
larvae  were  found  in  fast  water  along  the  remote  transect. 

Yet  as  the  summer  progressed,  the  number  of  animals  in  fast 
water  decreased,  whereas  abundance  near  the  shore  remained 
constant.  Thus,  the  relative  proportion  of  animals  in  slow 
water  near  shore  increased  through  the  summer.  This  micro- 
distributional  pattern  was  maintained  despite  the  much  greater 
nocturnal  drift  densities  observed  in  the  nearshore  area.  Hall 
et  al.  (I98O)  reported  a  similar  accumulation  of  larvae  of 
Tricorythodes  atratus  McD.  as  a  result  of  drift  into  pools 
from  riffle  areas  in  the  Mississippi  River.  Nevertheless, 
it  is  unclear  whether  the  increased  magnitude  of  drift  in  the 
shoreward  area  reflects  immigration  of  individuals  from 
faster  portions  of  the  river,  or  whether  it  represents  the 
desertion  of  the  slow-water  microhabitat  by  organisms  already 
there.  Bird  and  Hynes  (I98I)  have  also  observed  greater  drift, 
especially  by  Baetis  spp. ,  occurring  at  the  margin  of  a  stream 
than  in  the  centre.  They  attributed  this  to  animals  actively 
leaving  the  slow-water  areas.  Other  authors  have  also  reported 
increased  drift  in  response  to  slow  or  declining  water  levels 
(see  reviews  in  Armitage  1977  and  Walton  1980a). 
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In  autumn,  E.  inermis  larvae  were  concentrated  in  the  slow 
water  areas  near  shore.  Later  in  the  season,  however,  animals 
appeared  to  move  towards  the  river  centre.  The  adaptive  sign¬ 
ificance  of  this  seems  clear,  since  in  winter  ice  on  the  Pembina 
River  reaches  a  thickness  of  60  cm,  and  only  centre  portions 
of  the  river  maintain  a  flow  of  water.  Ulf strand  (1968) 
postulated  a  similar  phenomenon  to  account  for  sudden  changes 
of  abundance  of  Ephemerella  aurivillii  Bgtn.  in  a  north  Swedish 
river . 

This  movement  of  larvae  is  probably  accomplished  by 
individuals  crawling  over  or  through  the  substrate  rather  than 
by  drift.  Few  species  drift  abundantly  during  autumn  (Waters 
1972,  Ciborowski  1979)  even  though  benthic  densities  of  many 
species  are  maximal  at  this  time  (Hynes  1970)*  Drift  densities 
of  E.  inermis  larvae  in  October  in  the  Pembina  River  are  much 
lower  than  at  the  end  of  May  (Chapter  4).  Minshall  and 
Minshall  (1977)  found  that  E.  inermis  larvae  preferred  areas 
of  low  current  velocity  in  Mink  Creek,  Idaho,  a  stream  which 
apparently  does  not  become  ice-covered  in  winter.  Their  data 
suggest  that  the  response  to  current  velocity  was  weaker  in 
March  than  in  February.  Rabeni  and  Minshall  (1977)  found  E. 
inermis  larvae  to  have  no  preference  for  any  current  velocity 
in  Mink  Creek.  They  did  not  report  the  time  of  year  at  which 
their  study  was  conducted. 

There  appeared  to  be  no  clear-cut  associations  between 
abundance  of  either  species  and  concentrations  of  organic 
material  or  benthic  densities  of  other  invertebrates,  except 
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within  regions  of  similar  current  velocity.  Neither  Minshall 
and  Minshall  (19 77)  nor  Rabeni  and  Minshall  (1977)  found  a 
direct  relationship  between  detritus  accumulations  and  benthic 
densities  of  either  E.  inermis  or  B.  tricaudatus  larvae.  Amount 
of  FPOM  was  the  poorest  predictive  independent  variable  for 
regressions  grouped  by  replicate.  Hamilton  (1979)  reported 
that  this  organic  fraction  made  up  the  greatest  volume  of  gut 
contents  of  both  E.  inermis  and  B.  tri caudatus  larvae  in  three 
Alberta  rivers.  If  these  species  rely  on  the  same  food  source 
in  the  Pembina  River,  then  either  FPOM  is  sufficiently 
abundant  that  it  is  not  a  factor  important  in  controlling  the 
microdistribution  of  these  animals  or  the  larvae  are  so  active 
in  the  substrate  that  their  distribution  at  the  time  of 
collection  does  not  reflect  their  feeding  site. 

Corkum  et  al .  (1977)  observed  that  B.  tricaudatus  larvae 
preferred  organic  to  inorganic  substrates  in  laboratory 
experiments.  This  may  account  for  the  positive  association 
of  both  species  with  the  CPOM  fraction.  Alternatively,  the 
associations  of  both  species  with  all  three  independent  var¬ 
iables  may  simply  indicate  that  more  habitable  areas  (interstices) 
are  present  within  the  substrate  at  some  sampling  points  than 
at  others.  Interstices  act  as  traps  for  detritus  and  as 
refuges  from  the  current  for  all  invertebrates  (Rabeni  and 
Minshall  1977).  The  substrate  of  the  Pembina  River  consists 
of  smooth  cobbles  embedded  in  coarse  sand.  Suitable  refuges 
may  therefore  be  at  a  premium.  Regardless  of  the  relationships 
between  the  two  mayfly  species  and  the  three  independent  variables, 
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the  measurement  of  either  detritus  concentrations  or  density 
of  other  benthos  provides  a  means  of  statistically  reducing 
intersample  variability. 

Both  E.  inermis  and  B.  tricaudatus  larvae  are  vagrant 
species  capable  of  rapid  colonization  of  newly  available  substrates 
(Gore  1979 »  Shaw  and  Minshall  1980) .  Thus,  their  response  to 
the  varying  physical  characteristics  of  a  river  or  to  internally 
controlled  changes  in  habitat  preference  is  to  be  expected. 

The  microdistribution  of  both  species  appears  to  depend  largely 
on  the  constancy  of  prevailing  physical  conditions  in  the 
Pembina  River.  In  early  summer,  during  periods  of  stable  flow, 
both  species  predominate  in  faster  water.  Although  B.  tricaudatus 
larvae  maintain  their  position  in  these  fast  water  areas 
throughout  their  development,  E.  inermis  larvae  may  drift  to 
slower  water  as  they  reach  maturity.  Newly  hatched  larvae  of 
E.  inermis  are  most  abundant  near  shore  but  travel  to  deeper 
and  faster  water  in  late  autumn.  Flood  can  produce  a  rapid 
relocation  of  individuals  of  both  species,  probably  by  inducing 
animals  to  drift.  This  results  in  a  more  even  distribution  of 
animals  across  the  river.  The  distribution  of  both  species 
within  regions  of  similar  current  velocity  at  any  given  time 
is  associated  with  areas  of  high  detrital  concentrations, 
especially  the  coarse  fractions.  The  microdistribution  of 
either  species  clearly  depends  upon  responses  of  animals  to 
the  interaction  of  a  number  of  biotic  and  abiotic  factors 
(Rabeni  and  Minshall  1977)  but  the  primary  controlling  factor 
appears  to  be  the  flow  regime  of  the  Pembina  River. 
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CHAPTER  2 


INFLUENCE  OF  CURRENT  VELOCITY,  DENSITY  AND  DETRITUS 
ON  DEPARTURE  FROM  THE  SUBSTRATE  OF  TWO 
MAYFLY  SPECIES  (EPHEMEROPTERA) 
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Introduction 

The  microdistribution  of  lotic  aquatic  invertebrates  is 
largely  controlled  by  current  velocity  (Phillipson  1956, 

Edington  1965,  Madsen  1969,  Minshall  and  Minshall  1977,  Rabeni 
and  Minshall  1977),  substrate  composition  (Ulfstrand  1967, 

Cummins  and  Lauff  1969.  Higler  1975*  Hildrew  and  Townsend 
1977*  Williams  1978,  Barton  1980)  and  food  availability 
(Egglishaw  1964,  Mackay  1969*  Hildrew  and  Townsend  1976, 

Peckarsky  1980a) .  Drift  provides  an  effective  means  of 
removing  animals  from  microhabitats  in  which  one  or  more  of 
these  factors  becomes  unfavourable  (Corkum  et  al .  1977*  Walton 
et  al.  1977*  Corkum  and  Clifford  1980) .  Current  velocity  is 
of  major  importance  because  it  controls  the  distribution  of 
substrate  and  food  as  well  as  exerting  a  direct  influence  on 
benthic  invertebrates. 

Other  workers  have  postulated  that  drift  of  aquatic 
invertebrates  also  reflects  intraspecific  competition  for 
space  (see  reviews  in  Waters  1972,  Hildebrand  1974).  High 
benthic  densities  of  animals  may  increase  the  frequency  of 
interactions  between  individuals.  This  can  stimulate  benthic 
activity  and  result  in  either  a  greater  probability  of  accidental 
dislodgement  (Elliott  1967a,  McLay  1968) ,  or  the  active  desertion 
of  the  microhabitat  (Walton  et  al .  1977*  Hildrew  and  Townsend 
1980) .  Density-dependent  drift  occurs  when  animals  are 
abundant  enough  to  exceed  the  carrying  capacity  of  the 
microhabitat.  In  some  rivers,  however,  physical  perturbations 
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(e.g.  floods)  and  predation  may  maintain  population  densities 
below  this  level,  with  the  result  that  drift  is  density- 
independent  at  most  times  (Bishop  and  Hynes  1969a) . 

This  study  was  undertaken  to  examine  the  relative 
importance  of  current  velocity,  benthic  density  and  abund¬ 
ance  of  detritus  on  departure  tendencies  of  E.  inermi s  and 
B.  tricaudatus  larvae.  More  specifically,  I  attempted  to 
determine  whether:  (1)  larvae  exhibit  density-dependent 
drift  at  benthic  densities  typical  of  the  Pembina  River; 

(2)  presence  of  an  organic  substrate  (detritus)  modifies 
the  amount  of  space  available  to  animals  in  the  substrate; 

(3)  these  influences  are  of  greater  importance  in  controlling 
departure  rates  of  organisms  than  is  current  velocity. 
Departure  of  animals  from  the  substrate  was  monitored  in 
laboratory  streams  at  various  benthic  densities  under 
conditions  in  which  the  inorganic  substrate  was  augmented 
with  controlled  amounts  of  detritus.  Additional  experiments 
were  performed  to  estimate  the  influence  of  current  velocity 
on  departure.  The  response  of  dead  animals  to  current 
velocity  was  also  examined  to  discern  to  what  extent 
departure  is  the  result  of  accidental  physical  dislodgement 
of  animals. 


Methods 


Laboratory  Streams 

Experiments  were  conducted  in  May  (spring)  and  October 
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(autumn)  1979-  They  were  performed  in  three  elliptical, 
recirculating  artificial  streams  (channel  width  11  cm,  water 
depth  13  cm;  Fig.  10).  Water  in  each  stream  was  propelled 
by  a  motor-driven  paddlewheel .  Current  velocity  was  reg¬ 
ulated  by  varying  the  size  of  the  pulleys  that  coupled  the 
paddlewheel  with  the  motor.  The  streams  were  housed  within 
a  controlled-environment  room,  the  temperature  being  maintained 
at  prevailing  Pembina  River  temperatures  (Chapter  1). 

Experimental  Design 

Two  series  of  experiments  were  conducted  in  each 
season.  The  first  series  tested  the  influence  of  benthic 
animal  density  and  amount  of  detritus  in  the  substrate  on 
departure  of  E.  inermis  and  B.  tricaudatus  larvae.  A 
3  (density)  by  4  ( detritus)  level  factorial  stratified  random 
design  was  used,  incorporating  3  replicates  of  each  treat¬ 
ment.  Levels  corresponded  to  Jx ,  IX  or  3X  maximum  recorded 
abundance  in  spring  1978  bottom  samples.  The  fourth  detritus 
category  was  no  detritus  present  (Table  4) . 

The  second  series  of  experiments  examined  the  effect 
of  current  velocity  on  departure  of  larvae  from  the  sub¬ 
strate.  Three  replicates  were  conducted  for  each  of  5  mean 
current  velocities,  representing  the  range  of  current  found 
in  the  Pembina  River  up  to  a  distance  of  12  m  from  shore 
during  non-flood  flow  conditions.  Trials  with  live  animals 
were  all  conducted  with  3X  levels  of  density  and  detritus. 
Treatments  were  the  laboratory  equivalents  of  17*0,  24.7* 
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Table  4.  Animals  and  detritus  added  to  laboratory  streams.  All  confined 
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31.2,  40.8,  and  51*4  cm  s  ^  mean  velocity. 

Laboratory  equivalents  to  mean  velocity  in  the  Pembina 
River  were  estimated  by  taking  a  series  of  field  current 
velocity  readings  with  a  Guriy  current  meter.  Mean  velocity 
was  determined  by  taking  a  measurement  at  0.6X  stream  depth. 
Eight  points  were  selected  at  increasing  distances  from 
shore.  Mean  current  velocity  and  velocity  5.2  cm  above 
the  substrate  (equivalent  to  0.6X  laboratory  stream  depth) 
were  measured  at  each  point.  The  relationship  between  the 
two  measurements  was  determined  by  linear  regression  (Fig.  11 

Preserved  animals  from  earlier  experiments  were  used 
to  determine  the  effect  of  current  velocity  on  loss  from 
the  substrate  of  dead  organisms.  In  these  trials,  80 
E.  inermis ,  80  B.  tricaudatus  (spring)  and  80  B.  tricaudatus 
(autumn)  larvae  were  used.  Baetis  tricaudatus  larvae  were 
distinguished  by  size  and  maturity.  Those  from  spring 
trials  had  well-developed  wing  pads  and  head  widths  >  0.7 
mm.  Larvae  from  autumn  experiments  had  poorly-developed 
wing  pads  and  head  widths  of  0.4  -  0.6  mm.  These  larvae 
had  been  preserved  in  Kahle's  fluid  after  use  in  original 
experiments.  They  were  maintained  in  dechlorinated  tap 
water  for  1  h  prior  to  use. 

Several  workers  have  indicated  that  invertebrate 
predators  influence  drift  of  mayfly  larvae  (Corkum  and 
Pointing  1979,  Corkum  and  Clifford  1980,  Peckarsky  1980b, 
Walton  1980b).  Since  predators  are  ubiquitous  in  the 
Pembina  River,  they  were  used  in  all  trials  except  those 
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Figure  11.  Relationship  between  current  velocity  5.2  cm 

above  substrate  (X)  and  mean  current  velocity  (Y) 
in  Pembina  River.  Regression  equation  takes  the 
form  Y  =  -2.89  +  1.409X,  R  =  0.93* 
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concerned,  with  dead  animals,  .to  more  closely  approximate  field 
conditions.  I sogenoides  elongatus  (Hagen)  is  the  dominant  predator 
in  the  Pembina  River  and  was  added  to  laboratory  streams  in 
densities  that  approximated  average  abundance  of  all  predators 
in  the  field  (Table  4). 

Larvae  of  E.  inermis  were  not  used  in  autumn  trials, 
as  their  small  size  and  cryptic  habits  made  them  difficult 
to  collect  in  sufficient  numbers.  Experimental  densities 
of  B.  tricaudatus  larvae  in  autumn  were  doubled  over  spring 
levels  to  reflect  their  greater  abundance  at  that  time  in 
the  field. 

Execution  of  Experiments 

Detritus  for  all  experiments  was  collected  from 
the  Pembina  River,  sorted  to  remove  all  animals  > 3  mm 
in  length,  air-dried  and  weighed.  Particles  ranged  in  size 
from  0.420  mm  to  3*3^  mm.  The  required  amount  of  detritus 
for  each  trial  was  reconstituted  in  distilled  water  for 
24  h  prior  to  use.  This  provided  only  minimal  time  for 
bacterial  and  fungal  growth  to  occur.  Animals  were  hand- 
collected  from  the  river  and  maintained  in  aerated  water 
in  enamelled  trays  until  required. 

Detritus  was  placed  in  two  adjacent  11  x  30  cm 
trays  set  in  a  dry  laboratory  stream  and  covered  with 
scrubbed  substrate  from  the  Pembina  River.  Predators  were 
then  introduced.  The  stream  was  filled  with  refrigerated 
dechlorinated  tap  water,  with  the  paddlehweel  running,  and 
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left  for  1  hour.  The  appropriate  number  of  mayfly  larvae 
was  then  poured  into  the  stream.  All  larvae  eventually 
settled  on  the  substrate  instead  of  on  bare  portions  of 
the  stream.  Two  hours  after  introduction  of  the  mayfly 
larvae,  a  drift  net  was  placed  in  each  stream.  A  removable 
catchment  jar  on  each  net  was  replaced  at  30-minute  intervals 
for  the  next  4  hours.  The  jars  were  not  subsequently  changed 
until  the  end  of  the  light  period  of  the  experiment 
(spring,  8  h;  autumn,  5  h) .  The  lights  were  then  turned 
off  and  the  catchment  jars  changed  at  30-minute  intervals 
for  the  duration  of  the  experiment  (spring,  8  h;  autumn  10  h) . 

At  the  conclusion  of  each  experiment,  trays  containing 
substrate,  detritus,  predators  and  all  remaining  mayflies 
were  removed.  The  streams  were  drained  and  dried  prior 
to  the  beginning  of  the  next  trial. 

Results 

The  proportion  of  animals  leaving  the  substrate  during 
light  and  dark,  respectively,  was  determined  by  dividing 
the  number  of  animals  caught  in  drift  nets  during  relevent 
portions  of  each  experiment  by  the  number  present  at  the 
beginning  of  that  period.  Installment  of  drift  nets  produced 
a  noticeable  amount  of  disturbance  and  resulted  in  high 
initial  departure  rates.  For  this  reason,  the  first  3  h  of 
each  trial  were  not  considered  in  the  analyses.  All  data 
were  transformed  by  taking  the  arcsine  of  square  roots  of 
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proportions . 

The  relationship  between  proportion  of  animals  leaving 
the  substrate  and  benthic  density  was  tested  by  linear 
regression  analysis.  Separate  regression  lines  were  fitted 
for  each  level  of  detritus.  Analysis  of  covariance  was 
then  used  to  test  for  homogeneity  of  the  regression  coefficients. 
Proportions  of  animals  departing  under  different  detritus 
treatments  were  compared  by  1-way  analysis  of  variance. 

Curvilinear  regression  analysis  was  used  to  estimate 
the  relationship  between  proportion  of  animals  departing 
and  current  velocity.  Increasingly  higher-order  polynomials 
were  fitted  to  the  data  until  no  significant  improvement 
in  fit  (p  >  0.05)  was  observed. 

High  initial  departure  rates  of  both  live  and  dead 
insects  were  characteristic  of  all  trials.  In  all  cases, 
departure  stabilized  within  3  h  and  remained  fairly  constant 
for  the  duration  of  the  light  period  (Fig.  12) .  With  the 
onset  of  darkness,  departure  of  live  animals  increased 
substantially.  B.  tricaudatus  larvae  showed  a  more  pronounced 
nocturnal  response  in  spring  than  in  autumn  trials.  Departure 
of  live  animals  was  consistently  greater  than  erosion  of 
dead  insects.  The  total  proportions  of  B.  tricaudatus 
larvae  leaving  the  substrate,  exclusive  of  the  initial 
portion  of  experiments,  were  0.77  and  0.61  in  spring  and 
autumn,  respectively.  By  contrast,  only  0.08  and  0.07  of 
the  dead  animals  were  washed  away.  The  proportions  of  E. 
inermis  larvae  leaving  the  substrate  were  0.88  (live)  and 
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0 . 06  ( dead) . 

Of  the  24  regressions  performed  to  test  for  density 
dependence  (Table  5) »  only  one  produced  a  regression  coefficient 
significantly  different  from  zero  (p<0.05)»  I  interpret 
this  single  value  as  a  type  I  error.  Thus,  the  proportion 
of  animals  departing  from  the  substrate  was  independent  of 
benthic  density  for  both  species  under  all  levels  of 
detritus  in  light  and  in  darkness  (Fig.  13  and  14)  .  Analysis 
of  covariance  indicated  that  there  was  no  significant 
heterogeneity  of  regression  coefficients  (p>0.05,  Table 
6),  therefore,  the  interaction  of  detritus  and  benthic 
density  did  not  influence  departure  rates.  There  was  no 
significant  influence  of  amount  of  detritus  on  the  tendency 
for  animals  to  leave  the  substrate  (p>0.05,  Table  7;  Fig.  15)* 

In  spring,  departure  of  both  B.  tricaudatus  and  E. 
inermis  larvae  was  a  positive  linear  function  of  current 
velocity  (Fig.  16)  but  only  during  light  portions  of  the 
experiments.  During  darkness,  the  proportion  of  B.  tricaudatus 
larvae  leaving  the  substrate  was  inversely  related  to 

current  velocity.  Departure  of  E.  inermi s  larvae  was 

_1 

minimal  at  a  current  velocity  of  25  cm  s  and  greater  at 
both  higher  and  lower  velocities.  A  third-order  polynomial 
provided  the  best  fit  to  data  for  E.  inermis  larvae  in 
darkness  (Table  8).  Departure  of  B.  tricaudatus  larvae 
in  autumn  was  independent  of  current  velocity. 

Dead  animals  of  all  three  types  were  eroded  from  the 
substrate  in  much  lower  proportions  than  were  their  live 


1  ’  I  * 

■ 

* 

1 

« 


Table  5*  Significance  of  relationship  between  benthic  density  and  arcsine 
square  root  of  proportions  of  animals  departing  from  substrate  at 
various  levels  of  detritus.  *  indicates  p<0.05* 
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Figure  13.  Proportion  of  E.  inermis  larvae  departing  from  substrate  as  a  function 
of  benthic  density  during  light  (left)  and  dark  (right)  portions  of 
experiments.  Open  points  represent  two  observations. 
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Figure  14.  Proportion  of  B.  tri caudatus  larvae  departing  from  substrate  as  a  function 

of  benthic  density  during  light  (upper)  and  dark  (lower)  portions  of 
spring  (left)  and  autumn  (right)  experiments.  Open  points  represent  two 
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Figure  15.  Mean  proportion  and  95  percent  G.I.  of  B. 

tricaudatus  (upper)  and  E.  inermis  (lower)  larvae 
departing  from  substrate  as  a  function  of  detritus 
concentration  in  spring  (left)  and  autumn  (right) 
experiments.  Solid  intervals  represent  light 
portion  of  trials,  broken  intervals  represent  dark 
portioni 


Prop.  Departing 


Live  Dead 


Figure  16.  Proportion  of  live  (left)  and  dead  (right) 

larvae  departing  from  substrate  as  a  function 
of  mean  current  velocity.  Top,  E.  inermis 
(spring);  centre,  B.  tricaudatus  (springT; 
bottom,  B.  tricaudatus  (autumn).  Open 
points  and  broken  lines  represent  light 
portion  of  experiments.  Solid  points  and 
lines  represent  dark  portion. 
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c ounterparts  at  all  current  velocities.  Positive  linear 
relationships  were  found  between  current  velocity  and  de¬ 
parture  of  larvae  of  both  species  from  spring  experiments, 
but  only  during  light  periods  (Table  8).  Departure  of 
dead  insects  during  darkness  was  independent  of  current 
velocity  in  all  cases. 


Discussion 

Density-dependent  drift  responses  occur  when  the 
number  of  animals  present  in  the  substrate  exceeds  the 
space  available  (Waters  1972,  Walton  et  al.  1977,  Walton 
1978 »  Gyselman  I98O) .  Interactions  between  individuals 
at  sub-saturation  densities  may  produce  short-distance 
displacement  of  animals  but  seldom  result  in  drift  (Corkum 
1978a,  Wiley  and  Kohler  I98I).  Although  laboratory  densitie 
of  animals  considerably  exceeded  benthic  abundance  of  animals 
in  the  field  (Chapter  1),  departure  of  E.  inermis  and 
B.  tricaudatus  larvae  was  density-independent.  Thus, 
departure  of  these  organisms  in  the  Pembina  River  is  unlikely 
to  be  the  result  of  competition  for  space. 

Increasing  substrate  complexity  has  been  shown  to 
reduce  aggressive  intraspecific  interactions  between  lotic 
invertebrates  (Glass  and  Bovbjerg  1989)*  I  bad  anticipated 
that  should  interactions  between  larvae  be  responsible  for 
increased  departure  rates  of  animals  at  high  densities, 
the  addition  of  detritus  might  modify  this  response. 
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However,  since  all  drift  responses  were  density-independent, 
the  lack  of  a  significant  density-detritus  interaction  effect 
is  to  be  expected. 

Aquatic  invertebrates  deprived  of  an  adequate  food 
supply  usually  exhibit  a  density-independent  increase  in 
departure  rate  (Hildebrand  197^>  Keller  1976,  Bohle  1978). 

The  lack  of  any  departure  response  to  detritus  concentration 
in  my  experiments  may  indicate  that  relatively  sterile 
detritus  does  not  constitute  an  attractive  source  of 
nutrients  to  E.  inermis  or  B.  tricaudatus  larvae.  As  such, 
departure  rates  in  my  experiments  should  be  regarded  as 
maximal  rather  than  typical  of  the  field  situation. 

Behavioural  drift  refers  to  those  animals  that  have 
become  suspended  in  the  water  column  as  a  result  of  some 
change  in  their  internal  state  (Waters  1972).  Behavioural 
drift  may  be  active,  the  result  of  animals  releasing  them¬ 
selves  from  the  substrate,  or  passive,  whereby  activity 
or  positioning  renders  individuals  especially  susceptible 
to  accidental  dislodgement .  The  magnitude  of  the  behavioural 
component  of  drift  should  be  evident  as  the  difference 
between  transport  of  live  animals,  which  are  capable  of 
resisting  the  eroding  influence  of  the  current,  and  dead 
animals,  which  sire  deposited  and  then  eroded  from  the 
substrate  as  if  they  were  inert  particles. 

The  positive  linear  relationship  between  current 
velocity  and  departure  of  both  live  and  dead  animals  during 
the  light  portion  of  experiments  suggests  that  transport  is 


,  '  cil 

, 


67 


passive  a"t  'this  "time.  Live  animals  may  have  been  more 
prone  to  transport  than  dead  ones  because  they  selected 
exposed  areas  of  substrate  rather  than  lodging  in  inter¬ 
stices.  Corkum  et  al .  (1977)  reported  that  B.  tricaudatus 
larvae  inhabit  all  surfaces  of  substrate  particles,  and 
are  not  restricted  to  the  relatively  sheltered  undersides. 
Movements  of  the  animals  over  the  substrate  may  also 
contribute  substantially  to  the  increased  frequency  of 
erosion  of  live  insects  (Butz  1973)  especially  at  high 
current  velocities. 

Nocturnal  departure  rates  of  live  animals  in  spring 
were  greatest  in  trials  of  low  rather  than  high  current 
velocity.  This  suggests  an  added  active  behavioural 
response  on  the  part  of  the  organisms  during  darkness. 

High  incidence  of  nocturnal  drift  from  areas  of  reduced 
current  has  been  commonly  observed  (see  Walton  1980a). 

Corkum  et  al .  (1977)  found  that  B.  tricaudatus  larvae 
were  much  more  prone  to  leave  substrates  in  slow  water 
than  in  fast  water  under  both  light  and  dark  conditions. 

The  departure  of  dead  animals  in  the  dark  was  independent 
of  current  velocity,  presumably  because  most  animals  that 
lodged  in  the  most  exposed  areas  in  high  current  velocity 
trials  were  quickly  washed  away  early  in  those  experiments. 

Departure  of  B.  tricaudatus  larvae  in  autumn  experiments 
was  reduced  in  comparison  to  spring  trials  and  was  not 
influenced  by  current  velocity.  This  may  reflect  a  lower 
overall  level  of  activity  by  these  animals  as  well  as  a 
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broader  tolerance  for  reduced  currents.  Field  data  (Chapter 
1)  indicated  that  B.  tricaudatus  larvae  were  strongly 
concentrated  in  fast-water  areas  in  spring,  whereas 
distribution  was  more  even  with  respect  to  current  velocity 
in  early  autumn. 

Clearly,  the  marked  difference  in  departure  tendencies 
between  live  and  dead  animals  of  both  species  reflects  a 
major  behavioural  component  of  their  drift.  However,  of 
the  factors  investigated  in  this  study,  only  current  velocity 
exerted  a  significant  modifying  influence.  Passive  behavioural 
departure  produces  an  important  contribution  to  the  drift 
during  daylight,  but  it  is  the  active  nocturnal  component 
that  is  indicative  of  the  field  distribution  of  the  larvae. 
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CHAPTER  3 


DOWNSTREAM  AND  LATERAL  TRANSPORT  OF  DRIFTING  LARVAE 
OF  TWO  MAYFLY  (EPHEMEROPTERA)  SPECIES 
IN  THE  PEMBINA  RIVER 
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Introduction 


All  lotic  aquatic  invertebrates  occasionally  become 
dislodged  from  the  substrate  and  transported  downstream. 

The  rate  at  which  they  return  to  the  stream  bottom  depends 
on  the  size  and  shape  of  the  animals  (McLay  1970,  McKone 
1975)  as  well  as  the  current  velocity  (Elliott  1971a). 
Behavioural  actions  result  in  rapid  settling  rates  for 
certain  species  (Elliott  1971a) ,  especially  those  frequently 
appearing  in  the  drift  (Ciborowski  and  Corkum  I98O) . 

Most  studies  of  settling  behaviour  have  been  concerned 
with  the  distance  animals  are  carried  downstream;  but 
lateral  transport  of  suspended  animals  may  also  be  important 
in  determining  their  subsequent  benthic  microdistribution. 
Ulf strand  (1968)  postulated  that  cross-river  movements  of 
animals  could  produce  rapid  changes  in  microdistribution  of 
benthic  populations.  Lehmkuhl  and  Anderson  (1972)  reported 
that  animals  prone  to  drift  are  carried  from  riffles  to  back 
waters  at  the  stream  margin  during  floods.  However,  there 
is  no  experimental  work  that  evaluates  the  degree  of  lateral 
movement  exhibited  by  animals  released  into  the  water  column 

I  conducted  experiments  to  determine  settling  rates, 
downstream  movement  and  lateral  transport  of  the  mayflies 
Ephemerella  inermis  and  Baetis  tricaudatus  in  the  Pembina 
River.  To  determine  the  extent  to  which  the  animals  are 
capable  of  controlling  their  return  to  the  substrate,  I 
worked  with  both  live  and  dead  animals. 
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Methods 

Experiments  were  carried  out  during  the  last  2  weeks 
of  May  (spring)  and  the  first  2  weeks  of  October  (autumn) 
1979,  at  a  study  site  on  the  Pembina  River  (Chapter  1). 

The  site  was  mapped  with  respect  to  depth  and  mean  current 
velocity  prior  to  experiments  each  season.  A  base  line 
was  established  along  the  river  margin.  Eleven  transects, 

3  m  apart,  were  set  up  at  right  angles  to  the  base  line. 

I  recorded  mean  current  velocity  and  water  depth  at  1-m 
increments  along  each  transect  to  a  distance  of  at  least 
12  m  from  the  base  line.  Linear  regression  was  then  used 
to  estimate  the  relationship  between  distance  from  shore 
and  mean  current  velocity  for  each  season. 

Twelve  large  drift  nets  were  constructed  to  sample 
suspended  animals.  They  had  a  mouth  opening  of  1.0  x  0.5 
m  and  were  2  m  long.  Netting  was  composed  of  polyester 
sheer  drapery  material  with  a  pore  size  of  0.211  mm.  The 
nets  tapered  to  a  round  opening  (8  cm  diameter)  into  which 
were  sewn  jar  lids  with  centres  removed.  These  permitted 
rapid  attachment  and  removal  of  catchment  containers  (236 
ml  jars  with  bottoms  removed  and  replaced  with  netting) . 
Nets  were  anchored  by  means  of  aluminum  rods  driven  into 
the  substrate  to  sample  a  column  of  water  1  m  wide  and  up 
to  0.5  m  deep. 

The  twelve  nets  were  placed  contiguously  in  the  river , 
mouth  upstream,  in  a  row  extending  at  a  right  angle  to  the 
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base  line.  Settling  rates  were  determined  for  animals  at 
each  of  5  lines  established  parallel  to  the  base  line 
upstream  of  the  nets,  at  distances  of  2,  4,  6,  8,  and  10  m. 
from  shore.  Five  wooden  stakes  were  driven  into  the  sub¬ 
strate  along  each  line  to  serve  as  introduction  points 
for  the  insects  (Frontispiece).  Positioning  of  the  stakes 
varied  with  the  line  (Table  9)* 

Insects  were  hand-collected  from  the  river  and 
maintained  in  finger  bowls  containing  river  water  until 
required.  Those  designated  for  dead  animal  trials  were 
preserved  in  Kahle's  fluid  to  which  was  added  lignin  pink 
(a  vitalline  dye) .  These  insects  were  transferred  to  finger 
bowls  containing  river  water  30  minutes  before  they  were  used. 

A  trial  consisted  of  releasing  50  live  and  50  dead 
animals  of  each  species  into  the  water  just  above  the  stream- 
bed  at  a  single  stake.  Animals  were  released  by  placing 
them  in  a  118  ml  jar  containing  river  water,  immersing  the 
jar,  and  removing  the  lid.  Nets  were  submersed  for  two 
minutes  from  the  time  of  introduction  of  the  animals.  This 
was  judged  adequate  to  ensure  capture  of  suspended  animals 
at  the  slowest  current,  yet  minimized  inclusion  of  suspended 
detritus  and  extraneous  drifting  animals.  The  contents  of 
each  net  were  placed  in  separate  enamelled  trays  and  all 
appropriate  animals  were  removed  by  hand  and  counted. 

Only  large  animals  were  used  in  experiments  (head  widths 
>0.7  mm  in  in  spring,  0.4-  0.6  mm  in  autumn,  as  determined 
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by  later  measurements  in  the  laboratory) .  Ephemerella 
inermis  larvae  were  not  used  in  autumn  experiments  because 
their  small  size  and  cryptic  habits  made  them  difficult 
to  collect  at  this  time. 

To  estimate  the  number  of  extraneous  animals  caught 
during  experiments,  the  nets  were  used  to  sample  drift 
at  1800,  2400,  0900  and  1200  h  on  May  28-29  and  October  7-8 
1979*  Nets  were  submersed  for  from  2-5  minutes  in  spring 
and  from  15-25  minutes  in  autumn.  The  shorter  submersion 
period  of  nets  in  deepest  regions  was  used  to  ensure  that 
they  would  not  become  clogged  with  suspended  detritus. 

Net  contents  were  preserved  in  Kahle’s  fluid  and 
sorted  in  the  laboratory.  Head  widths  of  E.  inermis  and 
B.  tricaudatus  larvae  were  determined  with  an  ocular  micro¬ 
meter  under  a  dissecting  microscope  at  40X  magnification 
to  the  nearest  0.024  mm.  The  estimated  mean  number  of 
appropriately  sized  animals  caught  per  2  minutes  in  day¬ 
time  drift  samples  in  each  net  was  used  as  a  correction 
factor  for  settling  trials. 

Results 

'  Current  Velocity  Profile 

Mean  current  velocity  was  a  linear  function  of  distance 
from  the  river  margin  during  both  seasons  (Fig.  17  and  18). 
In  autumn,  however,  mean  current  velocity  appeared  to  become 
constant  at  distances  from  shore  greater  than  8  m.  Analysis 
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Current  Velocity  (cm/s) 
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Figure  1 7.  Relationship  between  mean  current  velocity 

and  distance  from  shore  in  May  1979-  Regression 
equation  takes  the  form  v  =  5-54  +  4.151y* 


.  » 


Current  Velocity  (cm/s) 
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Distance  From  Shore  (m) 


Figure  18.  Relationship  between  mean  current  velocity  and 

distance  from  shore  in  October  1979*  Regression 
equation  takes  the  form  v  =  -0.56  +  4.77 3y  for 
y<8  m,  and  v  =  36.3  for  y>8  m. 


of  variance  of  data  points  eight  and  more  meters  from 
shore  indicated  no  significant  increase  in  current  velocity 
occurred  past  this  point  (p>0.05).  Linear  regression  of 
all  data  points  in  spring  and  of  data  for  the  8  m  nearest 
shore  in  autumn  removed  O.77  and  0.84  of  the  variation, 
respectively. 

Relationship  Between  Settling  Rate  and  Mean  Current  Velocity 
Instantaneous  settling  rates  of  live  and  dead  animals 
of  each  species  for  each  line  were  determined  "by  the 
method  of  McLay  (1970).  The  relationship  between  number 
of  animals  suspended  in  the  water  column  (n  )  and  down- 

X 

stream  distance  (x)  between  point  of  release  and  the  nets 
is 

nx  =  n0e'rX  (1> 

where  n^  is  the  number  of  animals  released  and  r  is  the 
instantaneous  settling  rate.  A  value  of  r  was  estimated 
by  linear  regression  of  Ln(n  +1),  the  natural  logarithm 
of  the  number  of  animals  caught  in  all  nets  (corrected 
for  extraneous  animals,  and  adjusted  by  one  to  account 
for  zero  values),  against  x  for  each  transect. 

Mean  current  velocity  (v  )  for  each  line  (y  m 

1/ 

from  shore)  was  estimated  from  Fig.  17  and  18.  The  rel¬ 
ationship  between  r  and  v  was  then  estimated  by  regression. 

V 

Logarithmic  transformation  of  both  variables  best  described 
the  relationship  (Fig.  19) .  giving  regression  equations  of 
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the  form 

r  =  Glvy'B2  (2) 

where  13^  and  f3?  are  regression  constants.  Elliott  (1971a) 
and  Ciborowski  and  Corkum  ( 1 9 8 0 )  also  found  this  relationship 
to  best  fit  their  data. 

Settling  Capabilities 

To  compare  the  overall  settling  capabilities  of  animals, 
settling  capacity  (SC)  and  behavioural  settling  index  (Bl) 
were  calculated  (Ciborowski  and  Corkum  1980) .  Settling 
capacity  is  determined  by  integrating  equation  2  over 
the  range  of  current  velocities  tested  (approximately  10-50 
cm  s  ,  Fig.  19) •  The  relative  contribution  of  behaviour 
to  settling  ability  of  animals  (Bl)  is  estimated  by  taking 
the  complement  of  the  ratio  of  settling  capacities  of  dead 
to  live  animals . 

In  all  series  of  trials,  settling  rate  was  inversely 
related  to  current  velocity,  and  live  animals  exhibited 
greater  settling  rates  than  dead  ones  (Fig.  I9) .  In  spring 
experiments,  settling  rates  of  B.  tricaudatus  and  E.  inermis 
were  comparable  at  low  velocities,  but  E.  inermi s  larvae 
had  much  lower  settling  rates  at  higher  current 
velocities  than  did  B.  tricaudatus  larvae.  Large  B. 
tricaudatus  (spring  experiments)  had  a  greater  settling 
capacity  than  either  smaller  B.  tri caudatus  larvae  (autumn 
experiments)  or  E.  inermi s  larvae  (Table  10).  Behaviour  was 
of  more  importance  to  settling  of  small  than  large  B. 
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Table  10.  Regression  coefficients  of  relationship  between  instantaneou 
settling  rate  and  current  velocity,  and  measures  of  settling 
ability  of  animals. 
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tricaudatus  larvae;  however  both  had  a  value  of  BI  greater 
than  that  of  E.  inermis  larvae. 


Lateral  Transport  of  Suspended  Animals 


Diffusion  theory  predicts  that  suspended  particles 
should  disperse  normally  and  in  proportion  to  the  square 
root  of  the  distance  from  the  source  of  the  particles. 
Turbulence,  which  is  related  to  current  velocity,  should 
contribute  an  additional  diffusion  component  (Smith  1975) • 

I  used  multiple  regression  analysis  to  determine  how  the 
transverse  distribution  of  suspended  larvae  varied  with 
distance  from  point  of  release  (x)  and  with  current  velocity 
(v  ).  The  dependent  variable  was  the  standard  deviation  (s) 

y 

of  the  mean  distance  from  shore  of  suspended  animals  caught 

in  nets  for  each  trial.  Independent  variables  were  ^ 

and  v  .  Values  of  v  were  estimated  from  Fig.  17  and  18. 

V  V 

These  regression  equations  took  the  form 

0.5 


s  =  13^  +  B^vy 


(3) 


where  13^ ,  13^,  and  13^  are  regression  coefficients. 

Regressions  accounted  for  between  0.49  and  0.79  of 
the  total  variation  in  the  dependent  variable  (Table  11). 

All  13^  regression  coefficients  were  significantly  greater 
than  zero.  The  13^  coefficients  were  greater  than  zero  for 
all  analyses  except  that  for  E.  inermi s  (dead)  larvae  (p<  0.05). 
The  f3^  coefficient  consistently  explained  more  of  the  variation 
in  s  than  did  13^ .  Thus,  the  distance  animals  were  transported 
was  a  more  important  predictive  variable  than  was  mean  current 
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Table  11 . 


Regression  coefficients  of  relationship 
between  standard  deviation  of  mean  distance 
from  shore  of  suspended  animals,  distance 
downstream  from  point  of  release  (13,),  and 
current  velocity  (i3^).  13  represents  inter¬ 

cept.  Asterisks  inaicate^that  regression 
coefficient  is  significantly  greater  than 
zero;  *  p<.05,  **  p< . 01 ,  ***  p<.00l. 


% 

8  ^ 

G5 

R2 

E . 

inermis  (May) 

Live 

Dead 

-0.714 

-0.038 

0.570** 

0.243** 

0.0218* 

0.0028 

.56 

•59 

B. 

tricaudatus  (May) 

Live 

Dead 

-0.499 

0.151 

0.589*** 

O.O96*** 

0.0110* 

0.0028* 

•79 

•  56 

B. 

tricaudatus  (Oct.) 

Live 

Dead 

-0.494 

-0.243 

0.388* 

0.243*** 

0 . 0112* 
0.0061* 

.49 

.60 

* 

. 

♦ 
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velocity. 

Examination  of  the  regression  coefficients  in  Table  11 
indicates  that  dead  animals  should  be  less  dispersed  than  live 
animals  for  any  given  current  velocity  and  distance  below  a 
release  point.  However,  because  dead  animals  settle  over 
greater  distances  than  live  ones,  this  is  an  inappropriate 
comparison.  To  provide  a  better  comparison  of  potential 
dispersal,  both  between  live  and  dead  animals  as  well  as 
between  species,  I  determined  the  distance  required  for  one- 
half  of  each  sample  population  to  settle  (x  ^q) •  This 
can  be  calculated  from  the  formula  given  by  McLay  (1970), 

x  =  Ln(  1/0 . 50)r~'1' 

=  0.6931r-1  (4) . 


I  determined  values  of  x  ^  for  each  species,  using  values 
of  r  corresponding  to  hypothetical  releases  of  animals 
made  2  and  10  m  from  shore.  In  each  case,  r  was  calculated 
from  equation  2,  using  values  of  v  estimated  from  Fig.  17 
(Table  12).  The  estimates  of  x  represent  how  far  down¬ 
stream  of  the  release  point  animals  would  be  at  the  time 
one-half  had  settled.  The  mean  distance  from  shore  (y)  of 
the  25  animals  still  suspended  will  be  the  same  as  it  was 
when  they  were  initially  released,  but  they  will  now  be 
normally  distributed  with  a  standard  deviation  given  by 
equation  3*  I  constructed  frequency  distributions  for 
the  expected  transverse  distribution  of  suspended  animals 

by  substituting  x  <n  and  v  into  equation  3  to  obtain  expected 
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standard  deviations  (Table  12)  and  fitted  a  normal  distribution 
about  y  (Fig.  20) . 

Dispersal  of  suspended  animals  was  much  greater  along 
the  10  m  line  +han  along  the  2  m  line;  but  there 
were  no  consistent  trends  between  degree  of  transverse 
dispersal  of  live  and  dead  organisms.  Live  E.  inermis 
larvae  became  more  widely  dispersed  than  live  B.  tricaudatus 
larvae,  but  dispersal  differences  between  large  and  small 
B.  tricaudatus  larvae  were  minimal. 

Discuss! on 

Both  E.  inermi s  and  B.  tricaudatus  larvae  have  a  strong 
propensity  to  leave  the  substrate  and  become  transported 
by  the  river  current  (Chapters  1  and  2).  But  both  have 
developed  behavioural  adaptations  that  minimize  the  distance 
they  are  carried  downstream.  Ephemerella  larvae  have  little 
control  over  their  orientation  during  downstream  transport. 
These  animals  assume  a  characteristic  spread-leg  attitude 
and  exhibit  periods  of  undirected  swimming  activity  alternating 
with  periods  of  passivity.  Return  to  the  substrate  is 
largely  dictated  by  their  coming  into  chance  contact  with 
the  substrate  and  being  able  to  maintain  their  position 
there  (Elliott  1971a,  McKone  1 975 »  Luedke  and  Brusven  19 76, 
personal  observations).  Baetis  larvae  are  strong  swimmers 
and  are  thus  better  able  to  control  their  orientation 
during  transport,  especially  at  higher  current  velocities. 
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In  some  cases  these  larvae  are  capable  of  making  a 
directed  return  to  the  substrate  (Elliott  1971a,  McKone 
1975) •  For  both  species,  a  rapid  thigmotactic  response 
to  contact  with  the  substrate  is  of  primary  importance 
in  removing  animals  from  the  water  column. 

The  non-directed  swimming  actions  of  the  larvae 
have  the  additional  effect  of  dispersing  them  transversely 
in  the  water  column.  This  is  apparent  by  comparing  the 
regression  coefficients  of  live  and  dead  animals  of  each 
species  (Table  11).  At  any  given  distance  downstream  of 
a  point  of  release,  live  animals  will  be  more  widely 
dispersed  than  will  dead  ones.  But  the  greater  settling 
rate  of  live  animals  tends  to  reduce  differences  in  the 
overall  degree  of  transverse  dispersal  (Fig.  20) . 

This  lateral  transport  may  be  more  important  than 
the  degree  of  downstream  transport  that  animals  undergo. 

In  a  river  the  size  of  the  Pembina  (ca.  65  m  wide),  long¬ 
itudinal  changes  in  flow  structure  occur  over  distances 
much  in  excess  of  the  downstream  distance  an  animal  is 
likely  to  be  transported.  Thus,  it  is  improbable  that 
an  animal  will  be  carried  from,  for  example,  a  riffle  into 
a  pool  section  of  the  river.  In  contrast,  except  at  the 
river  centre,  current  velocity  changes  rapidly  across  the 
river.  Consequently,  animals  will  be  as  likely  or  more 
likely  to  experience  a  change  in  the  primary  physical 
factor  controlling  microhabitat  (current  velocity)  as  the 
result  of.  lateral  transport  than  because  of  downstream 
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Figure  20.  Frequency  distribution  of  suspended  animals  released 

2  and  10  m  from  shore  at  the  point  .at  which  one- 
half  settled.  Left  vertical  axis  indicates  animals 
remaining  from  an  initial  release  of  50.  Right 
vertical  axis  represents  probability  of  a  suspended 
animal  being  a  given  distance  from  shore.  Top,  E. 
inermis  (May);  centre,  B.  tricaudatus  (May); 
bottom,  3.  tri caudatus  (Oct.).  Solid  lines  represent 
live  animals;  broken  lines  represent  dead  animals. 


movement . 


The  instantaneous  distribution  of  suspended  larvae 
is  not,  by  itself,  an  accurate  measure  of  the  eventual 
distribution  of  larvae  in  the  substrate.  As  individuals 
become  increasingly  widely  laterally  dispersed,  they  come 
under  the  influence  of  increasingly  different  current 
velocities.  Those  insects  carried  towards  shore  will 
have  greater  settling  rates  than  those  carried  away  from 
shore.  Thus,  although  the  instantaneous  distribution  of 
suspended  animals  may  be  normal,  those  individuals  carried 
to  slower  water  will  have  a  greater  instantaneous  rate  of 
return  to  the  substrate  than  animals  carried  to  faster 
water.  The  distribution  curve  of  insects  that  actually 
do  settle  at  any  distance  (x)  below  a  release  point  should 
thus  have  a  mean  distance  from  shore  less  than  that  of 
their  release  point  (y) :  the  curve  should  be  skewed  towards 
the  river  centre.  The  skewness  of  the  distribution 
curve  will  be  determined  by  how  quickly  mean  current  velocity 
increases  as  a  function  of  distance  from  shore.  For  any 
given  animal  entering  the  water  column,  there  is  a  greater 
probability  of  returning  to  the  substrate  closer  to  shore 
than  of  becoming  reestablished  farther  from  shore. 

If  all  animals  in  a  population  were  equally  likely 
to  depart  from  the  substrate,  regardless  of  their  position 
in  the  river,  one  would  expect  this  population  to  become 
concentrated  along  the  river  margins,  given  ample  time. 
However,  both  B.  tricaudatus  and  E.  inermi s  larvae  are 
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more  prone  to  departure  from  slow  water  than  from  fast 
water  (Chapter  2).  This  may  normally  counteract  the  net 
effect  of  lateral  transport  of  suspended  animals.  Indeed, 
in  the  Pembina  River,  larvae  of  both  species  predominate 
in  fast  water  away  from  shore  during  non-flood  flow  conditions 
(Chapter  1).  When  the  river  floods,  however,  animals  do 
become  more  numerous  along  the  margins.  Flooding  results 
in  both  higher  current  velocities  in  the  river  as  well  as 
tending  to  increase  the  rate  of  erosion  of  animals  from  the 
substrate.  Lehmkuhl  and  Anderson  (1972)  also  reported  that 
flooding  carried  B.  tricaudatus  larvae  to  stream  margins. 

The  microdistribution  of  populations  of  animals  that 
are  prone  to  drift  possibly  represents  the  dynamic  equil¬ 
ibrium  achieved  by  a  balance  between  departure  and  settling 
tendencies  of  the  animals.  Prediction  of  the  benthic  distribution 
of  larvae  from  the  drift  would  require  determination  of  the 
joint  probabilities  of  departure,  downstream  movement  and 
lateral  transport.  Consideration  of  the  influence  of 
these  variables  on  microdistribution  may  provide  important 
insights  into  the  population  dynamics  of  these  species. 
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CHAPTER  4 


INFLUENCE  OF  DRIFT  ON  MICRODISTRIBUTION  OF  TWO 
MAYFLY  SPECIES  (EPHEMER OPTERA ) 
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Introduction 

Drift  is  an  important  aspect  of  the  population  dynamics 
of  many  lotic  invertebrate  species.  It  serves  to  remove 
individuals  from  certain  unfavourable  microhabitats  (Corkum 
et  al .  1977 »  Walton  et  al.  1977»  Walton  1978 >  Peckarsky 
1980b). and  provides  an  abundant  source  of  animals  that 
colonize  unexploited  areas  downstream  (Ulf strand  et  al . 

1974,  Townsend  and  Hildrew  1976,  Williams  and  Hynes  1976). 
Many  factors  can  induce  animals  to  depart  from  the  substrate 
(see  Thesis  Introduction).  The  manner  in  which  drifting 
individuals  return  to  the  substrate  has  also  been  elaborated 
(Chapter  3) •  Several  investigators  have  observed  that 
drift  can  result  in  seasonal  or  short-term  displacement 
of  significant  proportions  of  populations  (Ulf strand  1968, 
Lehmkuhl  and  Anderson  1972,  Hall  et  al .  1980) ,  but  the 
influence  of  drift  on  day-to-day  variation  in  relative 
distribution  of  animals  has  not  been  examined. 

In  this  chapter,  a  model  is  developed  that  examines 
the  role  of  daily  drift  in  maintaining  the  relative  distrib¬ 
ution  of  Ephemerella  inermis  and  Baeti s  tricaudatus  larvae 
in  an  ideal  river  with  flow  characteristics  similar  to 
those  of  the  Pembina  River.  In  the  Pembina  River,  these 
animals  are  usually  concentrated  in  fast  water  (Chapter  1). 
Nocturnal  departure  of  B.  tricaudatus  larvae  is  more 
pronounced  from  slow  water  than  from  fast  water  areas. 
Departure  of  E.  inermis  larvae  is  minimal  at  intermediate 
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water  velocities  (Chapter  2).  However,  individuals  of 
both  populations  entering  the  drift  tend  to  be  transported 
laterally  to  slower  water  as  well  as  downstream  (Chapter  3)« 
The  model  developed  herein  attempts  to  establish  whether  the 
benthic  distribution  of  these  populations  can  be  accounted 
for  by  simultaneously  considering  departure  and  settling 
capabilities  of  these  animals.  The  role  of  behaviour  in 
controlling  distribution  is  also  evaluated  by  comparing 
predicted  changes  in  distribution  of  live  individuals  with 
those  of  dead  animals,  which  behave  as  inert  objects. 

Methods 

Benthic  Distribution  as  a  Result  of  Drift 

Aquatic  invertebrates  can  change  position  in  a  river 
either  by  crawling  over  or  through  the  substrate  or  by 
entering  the  water  column.  For  the  purposes  of  the  model, 

I  will  initially  assume  that  the  latter  mechanism  is  the 
only  one  that  results  in  a  net  redistribution  of  animals. 
Drift  is  a  composite  process  made  up  of  three  events; 
departure  from  the  substrate,  transport  in  the  water  column 
(drift  proper),  and  return  to  the  substrate.  It  is  difficult 
to  separate  the  actual  moment  of  return  to  the  substrate 
from  the  period  of  transport.  Therefore,  I  will  refer 
to  these  two  components  together  as  settling. 

Consider  an  ideal  river  composed  of  a  discrete  number 
of  parallel  homogeneous  transects  of  arbitrary  width, 
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running  the  length  of  the  river.  Each  transect  has  an 
associated  known  mean  current  velocity.  Current  velocity 
is  minimal  at  river  margins  and  maximal  at  the  centre.  A 
fixed  proportion  of  animals  are  likely  to  enter  the  drift 
at  every  point  along  a  given  transect  each  day  because 
departure  is  a  function  of  current  velocity,  and  current 
velocity  is  uniform  along  the  length  of  each  cransect. 
There  is  also  a  fixed  probability  of  animals  returning  to 
the  substrate  in  the  original  or  any  other  transect  in 
the  river,  because  settling  and  lateral  transport  are  also 
functions  of  current  velocity.  If  mean  current  velocity 
varies  predictably  across  the  river,  then  probabilities  of 
departure  and  settling  can  be  estimated  by  considering  the 
position  of  any  particular  transect  relative  to  the  river 
margins . 


Let  =  (b^jb^tb^, . b?)  be  a  vector  where  each  b9 

represents  the  benthic  density  of  animals  in  transect  i  on 

a  given  day  0  (symbol  conventions  are  described  in  List  of 

Symbols).  Because  each  transect  is  uniform  along  its  length, 

b?  should  be  the  same  for  any  point  along  transect  i. 

Let  c. .  be  the  probability  that  an  animal  presently 

in  transect  i  will  be  located  in  transect  j  on  the  next  day. 

Were  it  possible  to  estimate  c. .  for  each  transect  of  the 

j 

river,  one  could  construct  a  probability  matrix,  C,  such 


that 
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Assuming  that  c^.  remains  constant  through  time,  then  C 
is  the  transition  matrix  of  a  Markov  chain,  and  one  can 
expect  B  ,  the  vector  describing  benthic  density  of  animals 
among  transects  after  t  days  to  be  given  by 


B 


t 


B°Ct 


(5) 


where  is  the  t'th  power  of  transition  matrix  C. 

Matrix  C  represents  the  probabilities  of  animals  in  one 
transect  being  located  in  others  after  one  day.  The  model 
requires  the  assumption  that  the  benthic  distribution  of  an 
animal  on  any  given  day  depends  only  on  its  location  on  the 
prior  day,  and  not  on  its  location  at  any  time  previous  to 
this . 


Because  not  all  animals  depart  from  the  substrate 
every  day,  let  P(t)  be  a  diagonal  matrix  where  each  diagonal 
element  p(t)^  represents  the  proportion  of  animals  that 
depart  from  the  substrate  during  time  interval  t.  Thus, 

Q(t)  is  the  proportion  of  animals  that  will  remain  on  the 
substrate , 


Q(t)  =  I  -  P(t)  (6) 

where  I  is  the  identity  matrix. 

Let  K  be  a  square  probability  matrix  in  which  each 

element  k. .  represents  the  probability  that  an  animal 
^  <J 

in  transect  i  that  has  left  the  substrate  settles  in 


transect  j . 
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Then  the  transition  matrix,  C,  can  be  determined  by 

C  =  P(t)K  +  Q(t)  (7) . 

That  is,  the  probability  of  animals  being  located  in  any 
transect  j  after  time  t  is  the  probability  of  departure 
from  the  substrate  in  any  transect  i ,  times  the  probability 
of  transport  from  transect  i  to  transect  j ,  plus  the 
probability  that  an  animal  in  transect  j  remains  there. 

I  assume  that  any  animal  that  does  drift,  does  not  drift 
more  than  once  per  day. 

The  Model  River 

In  Chapters  1  and  3»  I  examined  the  benthic  distribution, 
drift  and  settling  capabilities  of  live  and  dead  E.  inermis 
and  B.  tricaudatus  larvae  in  the  Pembina  River.  This 
river  is  approximately  65  m  wide  during  normal  flow  conditions, 
but  my  data  have  been  collected  primarily  from  a  region 
within  12  m  of  the  south  riverbank.  Within  this  area, 
both  mean  current  velocity  and  depth  are  linear  functions 
of  distance  from  shore  (Chapter  3).  Instead  of  extrapolating 
this  information  to  provide  estimates  of  flow  pattern  over 
the  entire  width  of  the  Pembina  River,  I  elected  to  base 
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predictions  on  a  hypothetical  river  with  flow  character¬ 
istics  analogous  to  those  of  the  sampled  portions  of  the 
Pembina  River.  The  model  river  is  26  m  wide  and  is  arbi¬ 
trarily  partitioned  into  1-m  wide  transects.  Depth  and 
current  velocity  between  0  and  12  m  from  the  south  bank 
of  the  model  river  increase  in  a  manner  identical  to  that 
estimated  for  the  Pembina  River.  Depth  and  current  velocity 
increase  at  twice  this  rate  between  0  and  6  m  from  the  north 
bank  of  the  model  river.  The  central  eight  meters  of  the 
model  river  are  of  uniform  depth  and  current  velocity  and 
are  equivalent  to  values  12  m  from  the  south  river  margin. 
Figure  21  illustrates  flow  characteristics  and  depth  profile 
of  the  Pembina  and  model  rivers  in  May  and  October  1979. 

Departure  of  Animals  From  the  Substrate 

In  Chapter  2,  regression  analysis  was  used  to  estimate 
the  relationship  between  departure  rates  of  live  and  dead 
E.  inermis  and  B  tricaudatus  larvae  and  mean  current  velocity. 
Estimates  of  departure  from  each  transect  of  the  model  river, 
during  day  and  at  night,  were  determined  by  substituting  the 
mean  current  velocity  of  each  transect  into  these  regression 
equations.  Where  departure  was  a  polynomial  function  of 
current  velocity,  mean  departure  rates  were  estimated  by 
integration.  Because  the  duration  of  laboratory  experiments 
was  different  than  field  photoperiods,  departure  rates  were 
corrected  to  a  16:8  L:D  photoperiod  in  May  and  to  a  12:12 
L:D  photoperiod  in  October. 
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The  model  permits  animals  to  dr 3  f  t  only  once  during  a 
24  h  period.  Let  p(L)^  and  p(D)^  be  the  proportions  of 
animals  departing  from  transect  i  during  day  and  night, 
respectively.  The  proportion  of  animals  in  transect  i 
that  remain  on  the  substrate  is 

Q.(L) ii  =  1  -  p(L)ii 
q(D)ii  =  1  -  p(D)ii 

The  proportions  that  depart  or  remain  during  24  hours  are 
given  by 

p(  2^0  ii  =  pCD).^  +  q(D)iip(L)ii  (8) 
q(24)ii  =  1  -  p(24)u  (9). 

Estimates  of  proportions  of  live  and  dead  animals  that 
depart  from  transects  of  the  model  river  are  presented 
in  Fig.  22 . 

Settling  of  Animals 

The  proportion  of  animals  that  settle  (f^)  over  a 
unit  length  of  any  transect  i  is  a  function  of  settling 
rate,  r,  and  downstream  distance,  x.  Let  m(x)^  be  the 
proportion  of  animals  suspended  in  the  water  column  over 
a  unit  length  in  transect  i,  x  meters  downstream  of  the 
point  of  entry  into  the  water  column.  When  x  =  0, 
m(x)^  =  The  proportion  of  animals  suspended  x 

meters  downstream  of  the  entry  point  is  given  by 
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Figure  22.  Proportion  of  animals  that  depart  from  substrate 

in  model  river  during  light  (solid  lines)  and 
darkness  (broken  lines).  A,  E.  inermis  (May); 

B,  B.  tricaudatus  (May);  G,  B.  tricaudatus 
(October!.  Figures  above  represent  departure 
of  live  individuals;  figures  below  represent 
departure  of  dead  animals . 
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m(x)^  =  m(0)^e~rix 


(10) 


(see  Chapter  3).  Therefore,  the  proportion  that  settle 
over  a  unit  of  downstream  distance  is  given  by 


f(x^)  =  m(x-l)^  -  m(x-l)^e  ri 

=  m( x-1 ) i ( 1  -  e_ri )  ( 11 )  . 


Settling  rate  is  itself  a  power  function  of  current 
velocity  (v)  such  that  for  any  transect  i, 


(12) 


where  f3^  and  f32  are  regression  constants  (see  Chapter  3)* 
However,  current  velocity  is  a  linear  function  of  distance 
from  shore  (y)  except  at  the  river  centre  (see  Fig.  21). 
Thus ,  for  any  y^ , 


(13) 


where  13^  and  13  2  are  regression  constants. 

The  lateral  distribution  of  suspended  animals  about 
their  point  of  introduction  into  the  water  column  (y^)  is 
normal  with  the  mean  given  by  distance  from  shore  (y^), 
and  standard  deviation  s^  (Chapter  3) •  The  standard  deviation 
is  dependent  on  both  downstream  distance  travelled  (x) 
and  current  velocity  (v)  such  that  for  animals  suspended 
in  any  transect  i , 


(14) 


■  • 
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where  and  are  regression  constants.  But  because 

current  velocity  is  a  linear  function  of  distance  from  shore, 

si  =  °3  +  B4X°'5  +  (15). 

Data  on  settling  rates  and  lateral  dispersal  of 
suspended  live  and  dead  E.  inermis  and  B.  tricaudatus 
larvae  (Chapter  3)  were  reanalysed  to  express  r  and  s 
as  functions  of  distance  from  shore  (Table  13).  An 
estimate  of  r^  and  s^  for  each  transect  i  of  the  12  m  from 
the  south  shore  of  the  model  river  was  determined  by  sub¬ 
stituting  the  distance  from  shore  of  the  midpoint  of  each 
transect  into  the  regression  equations.  Estimates  for 
each  transect  between  0  and  6  m  from  the  north  shore  were 
determined  by  substituting  twice  the  midpoint  distance  into 
the  regression  equations.  Values  for  the  eight  central 
transects  were  calculated  by  setting  y^  =  12  and  sub¬ 
stituting  into  the  equations  (Fig.  23) . 

Probabilities  of  transport  of  animals  from  one  transect 
to  all  others  were  generated  according  to  an  APL  computer 
program  (Appendix  1).  The  program  was  designed  to  determine 
proportions  of  animals  settling  at  increasing  1-m  steps 
downstream  of  the  given  point  of  entry ,  y^ ,  into  the  water 
column.  The  proportion  of  animals  not  settling  at  each 
step  was  partitioned  according  to  a  normal  probability 
distribution  generated  around  a  mean  (y^)  and  a  standard 
deviation  ( s. ) .  Proportions  returning  to  the  substrate  at 
each  step  were  summed  until  the  total  fraction  that,  had 


' 


, 


Table  13*  Regression  coefficients  and  coefficients  of  determination  of  functions 

estimating  instantaneous  settling  rate  and  measures  of  lateral  dispersal 
of  suspended  animals  in  model  river. 
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Figure  23*  Instantaneous  settling  rate  (  concave^  curves  , 

left  axes)  and  measures  of  lateral  dispersal 
(convex  curves,  right  axes)  of  live  (solid 
lines)  and  dead  (broken  lines)  larvae  suspended 
in  model  river.  A,  E.  inermis  (May);  B,  B. 
tricaudatus  (May) ;  C,  B.  tricaudatus  (October) . 


■ 


settled  exceeded  0.999*  The  remainder  was  added  to  the 
transect  from  which  animals  originated.  Summing  the  pro¬ 
portions  of  animals  suspended  at  the  end  of  each  step 


provided  an  estimate  (w.  .)  of  the  expected  contribution 

of  animals  originating  in  transect  i,  to  the  drift  at 

other  transects,  j  .  Each  run  of  the  program  produced  a 

vector  =  (k^k^k^. . k26^  ^h  each  kj  the 

probability  of  an  animal  originating  in  transect  i  settling 

in  transect  j,  and  a  vector  V\h  =  (w1,w2»w^, . w26^  ’ 

with  each  w.  the  sum  of  proportions  of  animals  originally 
J 

suspended  in  transect  i  that  became  part  of  the  drift  in 
transect  j .  Matrices  K  and  W  were  then  defined  as 


K26  /  \  “26  / 

Matrix  K  was  used  in  equation  7.  Matrix  fa  was  used  in 
determining  estimates  of  numbers  of  animals  that  would  be 
caught  in  drift  nets  placed  across  the  model  river  (see 
below) .  K  and  W  were  estimated  for  both  live  and  dead 
E.  inermis  and  B.  tricaudatus  larvae.  Matrices  are 


tabulated  in  Appendix  2. 


/ 
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Benthic  Distribution  of  Animals  in  the  Model  River 

Data  on  benthic  microdistribution  of  E.  inermis  and 
B.  tricaudatus  larvae  in  the  Pembina  River  were  collected 
throughout  the  ice-free  period  in  1979  (Chapter  1).  Five 
replicate  Hess  samples  were  taken  at  10-15  day  intervals 
at  distances  of  3  and  11  m  from  the  south  bank  until  mid- 
June,  and  at  3>  7»  and  11  m  from  the  south  bank  thereafter 
until  November. 

Benthic  densities  in  the  model  river  were  based  on 
field  benthic  densities  of  only  large  E.  inermis  and  B. 
tricaudatus  larvae  from  May  and  October  (Chapter  2). 

Numbers  of  large  animals  in  each  sample  within  series  of 
replicates  on  each  date  were  adjusted  by  regressing  benthic 
densities  against  numbers  of  all  other  animals  taken  in 
each  sample  (Chapter  1)  to  reduce  sample  variances  (Sokal 
and  Rohlf  19&9) • 

The  relationship  between  adjusted  benthic  density  of 
each  population  and  distance  from  the  south  river  margin 
on  each  date  was  determined  by  piecewise  linear  regression. 
Benthic  densities  of  animals  in  the  model  river  were  then 
estimated  from  these  regression  equations  (Table  14). 
Densities  between  13  and  20  m  from  the  south  shore  were 
the  same  as  at  the  12  m  point.  Benthic  densities  from 
0  to  6  m  from  the  north  bank  were  obtained  by  doubling 
the  slope  of  the  regression  equations. 


■ 
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Table  14.  Regression  coefficients  and  coefficients  of 

determination  of  functions  describing  benthic 
distribution  of  animals  in  model  river. 


Species 

Date 

Intercept 

Slope 

R2 

E. 

inermis 

May  15 
May  28 

-13.00 

3.25 

7.00 

8.25 

0 . 86 
0.80 

B. 

tricaudatus 

May  15 
May  28 

-13.49 

-23.53 

6.72 

8.50 

0.41 

O.63 

Oct .  4 

Oct  15 

-16 . 62 
-115.75 

20.01 

41 . 30 

0.63 

0.46 

* 
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Comparison  of  Observed  and  Predicted  Changes  in  Benthic 
Distribution 

All  settling  and  departure  estimates  were  based  on 
the  results  of  experiments  conducted  with  large  animals 
collected  in  May  or  October  1979*  Benthic  samples  that 
most  closely  corresponded  with  these  periods  were  taken 
on  May  1 5  and  28,  and  October  4  and  15,  1979. 

Predictions  of  relative  changes  in  benthic  density 
in  the  model  river  were  made  for  13  and  11  day  intervals, 
respectively,  to  correspond  to  the  time  elapsed  between 
samples  taken  in  the  Pembina  River.  Predictions  were 
made  using  equation  5*  For  May  trials,  the  initial 
distribution  of  animals  (B^)  was  the  vector  of  benthic 
densities  estimated  from  May  15  benthic  samples.  Initial 
distribution  (B^)  for  October  trials  was  that  estimated 
from  October  4  samples.  Dead  populations  were  given  the 
same  initial  distributions  as  their  live  counterparts. 

Comparison  of  Observed  and  Predicted  Drift  of  Animals 

Observed  Drift--  Drift  of  animals  in  the  Pembina 

River  was  studied  on  May  28-29  and  October  7-8,  1979 

(Chapter  3).  Large  drift  nets  filtered  all  animals  from 

12  1-m  wide  transects  extending  from  the  south  river  bank. 

Samples  were  taken  at  0900,  1200,  1800  and  2400  h.  All  E. 

inermis  and  B.  tricaudatus  larvae  were  measured  and  only 

large  individuals  were  used  in  subsequent  analyses.  Drift 

-3 

was  expressed  as  drift  density  (number  of  animals  m  water 
filtered) . 


■ 

' 
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Drift  of  E.  inermis  and  B.  tricaudatus  larvae  is 
aperiodic  during  daylight  hours  but  not  at  night  (see 
Fig.  8,  Chapter  1).  Maximum  drift  of  E.  inermi s  larvae 
occurs  near  midnight.  Baetis  tricaudatus  larvae  exhibit 
at  least  two  nocturnal  peaks  of  drift,  the  major  one 
occurring  just  prior  to  dawn.  Average  daytime  drift 
density  of  each  population  in  each  transect  was  determined 
by  dividing  the  total  number  of  animals  caught  in  the  three 
daytime  sampling  periods  by  the  total  volume  of  water 
filtered  by  each  net.  Average  nocturnal  drift  was  estimated 
by  dividing  the  number  of  animals  caught  in  each  net  at 
2400  h  by  the  volume  of  water  filtered  at  that  time  and 
multiplying  these  values  by  the  difference  between  drift 
at  2400  h  and  average  nocturnal  drift  in  Fig.  8. 

Predicted  Drift--  The  total  numbers  of  animals 

2 

expected  to  depart  from  one  m  of  substrate  in  transect  i 
of  the  model  river  during  day  (L)  and  at  night  (D)  are 
given  by 

n(L)i  =  bip(L)ii  (16) 

=  bip(D)ii  (17) 

where  b^  is  the  benthic  density  of  animals  in  transect 
i  and  p^  is  the  proportion  that  depart  from  the  substrate. 
Since  matrix  W  represents  the  sum  of  proportions  of  animals 
departing  from  transect  i  that  would  be  carried  over  a 
point  in  transect  j  ,  the  total  number  of  animals  that 
should  be  caught  in  drift  nets  placed  across  the  transects 


during  day  (N  (L))  and  night  (N  (D))  would  be 

iiid-A.  rna.x 

Nmax(L)  =  N(L)W  (18) 

Nraax(D)  =  N(D)M  (19). 


However,  the  concentration  of  animals  in  the  water  above 
2 

each  m  of  each  transect  also  depends  on  the  depth  of 
water  (z^)  in  each  transect  i.  Thus,  the  concentration 
of  drifting  animals  in  each  transect  is  given  by 


n' 

n' 


max 


(L) 


max 


(D) 


i 

i 


n  (L)  .  /z. 
max  v  i 

n  (L ) .  /z. 
max'  i'  i 


(20) 
(21)  . 


Volumes  of  water  passing  points  in  each  transect 
of  the  model  river  were  estimated  using  current  velocity 
and  depth  estimates  from  Fig.  21  and  assuming  photoperiods 
of  1 6  s  8  L:D  in  May  and  12:12  L:D  in  October.  Predicted 
drift  density  during  day  and  at  night  was  determined  by 
dividing  the  total  diluted  number  of  animals  (N'  ) 

predicted  to  drift  past  a  point  during  day  and  at  night 
by  the  estimated  volume  of  water  passing  each  point. 
Expected  drift  density  of  live  and  dead  B.  tricaudatus 
larvae  was  determined  using  benthic  distributions  estimated 
on  May  28  and  October  4.  Expected  drift  of  live  and  dead 
E.  inermis  larvae  was  determined  by  using  estimated 


benthic  distributions  on  May  28. 


. 


■ 

' 
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Mean  Daily  Downstream  Distance  Travelled  by  Animals 

The  sum  of  proportions  of  animals  that  drift  past 
a  point  in  a  river  is  the  product  of  the  proportion  departing 
from  the  substrate  and  the  mean  distance  travelled  (Elliott 
1971a).  Thus, 

w  —  m ( 0 ) x  (22) . 

When  m(0)  =  1,  as  assumed  in  the  settling  computer 
program,  this  equation  reduces  to 

x  ~  w  (23) . 


Since  W  is  the  matrix  of  the  sums  of  proportions  of 
suspended  animals  that  are  transported  from  any  transect 
i  to  each  transect  j ,  the  mean  downstream  distance  travelled 
by  animals  that  depart  from  any  transect  i  (x^)  is  the 
i '  th  row  total  of  W,  i.e. 


x.  =  w.  1  +  w.  0  +  w.  ~  + 
i  ll  i2  i3 


+  wi26 


(24) 


and 


X  =  (x1,x2,x^, 


x26). 


The  mean  distance  travelled  by  all  animals  in  the  population 
that  drift  during  one  day  (x^)  is 


XN  N ( 24) I 


N(24)X 


(25) 


where  N(24)  is  the  vector  of  the  number  of  animals  that 
depart  from  the  substrate  over  24  hours  and  I  is  the 
identity  vector. 


Ill 


The  mean  downstream  displacement  of  the  entire  population 
as  a  result  of  drift  is 


N(24)X 

B°I 


(26) 


where  B°  is  the  vector  of  benthic  density  of  animals 

across  the  model  river  on  day  0.  It  is  important  to 

note,  however,  that  because  n(24).  =  b.p(24)..,  both 

1  i  n 

xN  and  xB  are  dependent  on  the  benthic  distribution  of 
animals . 

Downstream  distances  travelled  by  drifting  animals 
and  displacement  of  populations  were  calculated  for  live 
and  dead  E.  inermis  and  B.  tricaudatus  larvae,  based  on 
estimates  of  observed  benthic  distributions  on  May  28  and 
October  4  1979* 


Results 

Benthic  Distribution 

The  observed  benthic  distribution  of  animals  changed 
minimally  over  the  time  intervals  examined  in  May  and 
October.  In  contrast,  the  model  predicted  large-scale 
transport  of  individuals  of  all  populations  towards  shore 
(Fig.  24, 25 >26).  In  all  cases,  predicted  densities  of 
animals  fell  outside  the  95  percent  confidence  intervals 
of  observed  densities  of  organisms.  The  model  also  indicated 
that  densities  of  animals  should  become  greater  adjacent  to 
the  north  bank  than  the  south  bank  of  the  model  river. 
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Ephemerella 


A 


Figure  24.  Observed  (solid  lines)  and  predicted  (broken) 

lines)  distribution  of  E.  inermis  larvae  in 
model  river.  Vertical  lines  represent  95  percent 
C.I.  of  observed  densities  of  larvae.  A,  May  15 
1979;  B,  May  28  1979  (live  animals);  C,  May  28 
1979  (dead  animals). 


Benthic  Density,  no. 
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Baetis 


May  28 


C 


Figure  25.  Observed  (solid  lines)  and  predicted  (broken 

lines)  distribution  of  B.  tricaudatus  larvae 
in  model  river.  Remaining  explanation  as  in 
Fig.  24. 
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Figure  26.  Observed  and  predicted  distribution  of  a 

tricaudatus  larvae  in  model  river.  A, 

October  4  1979;  B,  October  15  1979  (live 
animals);  C,  October  15  1979  (dead  animals). 
Remaining  explanation  as  in  Fig.  24. 
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Changes  in  distribution  predicted  for  live  animals 
were  more  pronounced  than  for  dead  populations  of  the 
same  species.  A  greater  degree  of  lateral  redistribution 
was  evident  in  May  than  in  October  populations  of  B. 
tricaudatus  larvae.  In  May,  changes  in  benthic  distribution 
of  the  E.  inermi s  population  were  more  noticeable  than 
changes  in  distribution  of  the  B.  tricaudatus  population. 

Drift  Density 

Predicted  drift  densities  of  populations  of  dead 
animals  were  generally  much  lower  than  predicted  drift 
densities  of  corresponding  live  populations. 

Observed  drift  densities  of  E.  inermi s  larvae  were 
much  greater  in  May  than  in  October  (Fig.  27) .  In  May, 
predicted  drift  densities  of  live  animals  corresponded 
closely  to  values  observed  in  the  Pembina  River  both  during 
day  and  at  night. 

Observed  nocturnal  drift  densities  of  B.  tricaudatus 
larvae  were  slightly  greater  in  May  than  in  October  (Fig. 
28).  The  reverse  was  true  for  daytime  drift  density  values. 
The  model  markedly  underestimated  May  drift  densities  of 
live  B.  tricaudatus  larvae.  October  predictions  of  drift 
density  were  too  high  during  daytime  and  too  low  at  night. 
All  predictions  of  drift  densities  of  dead  animals  were 
underestimates  of  drift  density  of  live  animals  in  the 


Pembina  River. 


V 

. 

' 


11<: 


B 


=t — r^'l 

14  18  22  26 


DISTANCE  FROM  SHORE,  M 

Figure  27.  Observed  (solid  lines)  and  predicted  (broken  lines) 

drift  density  of  E.  inermis  larvae  in  Pembina  and 
model  rivers,  respectively,  in  May  (left)  and  Oct¬ 
ober  (right).  A,  daytime  drift  density  (live 
animals  above,  dead  animals  below);  B,  night¬ 
time  drift  density  (live  animals  above,  dead 
animals  below) .  Note  that  scales  of  vertical 
axes  differ. 
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Figure  28.  Observed  and  predicted  drift  density  of  B. 

tri caudate s  larvae  in  Pembina  and  model 

.  Remaining  explanation  as  in  Fig.  27. 


rivers 


118 


Distance  Drifted  Downstream 

Dead  individuals  of  both  species  that  departed  from 
the  substrate  were  transported  at  least  twice  as  far  down¬ 
stream  as  their  live  counterparts  (Table  15).  However, 
because  live  animals  had  higher  departure  rates,  the 
daily  downstream  displacement  of  entire  populations  was 
greater  for  live  than  for  dead  animals.  Downstream  dis¬ 
placement  of  the  E.  inermis  population  was  over  twice  as 
great  as  for  the  B.  tricaudatus  population.  Displacement 
of  B.  tricaudatus  populations  (live  animals)  were  almost 
identical  in  May  and  October. 

Discussion 

Results  of  the  model  runs  indicate  that  drift  is 
responsible  for  transport  of  large  numbers  of  animals  from 
the  river  centre  to  near-shore  regions  as  well  as  downstream. 
This  effect  was  observed  for  both  populations,  and  was  most 
noticeable  in  live  animals.  Test  runs  of  the  model  using 
theoretical  data  (Appendix  3)  indicate  that  animals  will 
become  concentrated  at  stream  margins  whenever  lateral 
dispersal  values  at  shoreward  transects  are  close  to  zero, 
regardless  of  individuals'  departure  and  settling  capabilities. 
My  estimates  of  lateral  dispersal  for  live  animals  in  both 
populations  included  at  least  one  such  value.  Thus,  the 
model  suggests  that,  given  enough  time,  drift  should 
result  in  all  E.  inermis  and  B.  tricaudatus  larvae  eventually 
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Table  15- 


Predicted  mean  downstream  distance  travelled 
by  animals  that  drift  (x^)  and  net  downstream 
displacement  of  populations  (3L)  in  model 
river.  u 


Species 


Condition  x^ 


E.  inermis  (May) 


B.  tricaudatus  (May) 


Live 

15.10 

14.89 

Dead 

37.29 

5.59 

Live 

8.25 

6.44 

Dead 

17.70 

3.24 

Live 

8.99 

6.43 

Dead 

26.00 

1.76 

B.  tricaudatus  (Oct.) 


. 
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becoming  located  in  and  remaining  in  those  transects  closest 
to  river  margins.  Estimates  of  change  in  distribution  of 
dead  organisms  were  less  pronounced,  primarily  because 
dead  animals  were  less  prone  to  depart  from  the  substrate. 
Reduced  departure  rates  tend  to  slow  the  rate  of  change  in 
distribution,  even  when  settling  rates  of  suspended  animals 
are  low  (Appendix  3)« 

Differences  in  the  rate  of  change  in  distribution 
among  populations  of  live  animals  were  also  largely  due 
to  differences  in  daily  departure  rate.  It  appears  that 
behavioural  mechanisms  that  result  in  reduced  departure 
rates  from  fast-water  areas  are  inadequate  to  compensate 
for  the  increased  degree  of  lateral  transport  associated 
with  behavioural  mechanisms  that  facilitate  a  rapid  return 
to  the  substrate.  Indeed,  the  influence  of  behaviour 
results  in  more  pronounced  lateral  displacement  of  animals 
than  would  occur  if  individuals  departed  and  settled  in 
a  completely  passive  manner. 

Examination  of  predicted  downstream  displacement 
estimates  shows  a  similar  situation.  Even  though  dead 
animals  that  became  suspended  in  the  water  column  were 
transported  greater  distances  downstream  than  their  live 
counterparts ,  higher  departure  orates  of  live  animals 
resulted  in  greater  net  downstream  displacement  of 
populations . 

Because  there  was  such  a  large  discrepancy  between 
observed  and  predicted  changes  of  distribution,  it  is  tempting 
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to  suggest  that  the  results  are  an  artifact  of  my  estimates 
of  departure,  settling  or  lateral  transport.  Yet,  model 
estimates  of  the  magnitude  and  spacial  pattern  of  drift, 
which  were  based  on  the  same  data,  were  surprisingly  close 
to  drift  observed  in  the  Pembina  River.  Although  I  had 
expected  the  model  to  overestimate  amount  of  drift  because 
derived  departure  rates  of  live  animals  were  almost  certainly 
overestimated  (see  Chapter  2),  the  most  pronounced  deviations 
were  underestimates.  It  is  possible  that  overestimates  in 
departure  rate  were  compensated  for  by  equally  erroneous 
overestimates  of  settling  rate.  All  determinations  of 
settling  rate  were  made  during  daylight  hours  (Chapter  3) > 
whereas  most  drift  occurred  at  night,  especially  in  May. 
Luedke  and  Brusven  (I976)  and  Walton  (1978)  found  that 
drifting  animals  tended  to  have  lower  settling  rates  at 
night  than  during  the  day.  However,  their  experiments 
were  carried  out  in  very  shallow  laboratory  streams  (depth 
less  than  6  cm);  thus,  their  observations  may  not  be 
applicable  to  all  portions  of  a  river.  Elliott  (1971a) 
found  no  differences  between  daytime  and  nocturnal  settling 
rates  of  animals  released  into  the  water  column  of  a  small 
English  stream. 

An  alternative  explanation  is  that  some  animals 
engage  in  multiple  drift  episodes.  In  such  an  event,  the 
model  would  underestimate  drift  because  my  estimates  of 
departure  assumed  that  animals  drifted  not  more  than  once 
every  24  hours.  It  should  be  noted,  however,  that  any  true 
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underestimate  of  the  number  of  animals  transported  in  the 
water  column,  or  of  the  downstream  distance  they  are  carried, 
would  produce  a  corresponding  underestimate  of  the  amount 
of  lateral  redistribution  of  the  population.  For  this 
reason,  I  feel  that  the  predicted  changes  in  distribution 
are  representative  of  the  direction,  if  not  the  exact 
magnitude,  of  the  transportive  effects  of  drift  on  the 
distribution  of  E.  inermis  and  B.  tricaudatus  populations. 

Because  there  is  such  a  large  discrepancy  between 
observed  and  predicted  benthic  distribution  of  larvae, 
it  is  clear  that  drift  alone  cannot  account  for  the  micro¬ 
distribution  of  these  populations  in  the  Pembina  River. 

The  benthos  must  undergo  some  form  of  compensatory  benthic 
movement  that  counteracts  the  transportive  effects  of  drift. 
Although  there  is  no  experimental  evidence  that  aquatic 
invertebrates  are  capable  of  directed  lateral  movements 
in  rivers,  indirect  information  from  several  sources  lends 
credence  to  the  possibility.  Seasonal  onshore-offshore 
movements  of  mayfly  populations  in  lentic  habitats  have 
been  reported  by  Gibbs  (1973>  1979) »  Corkum  ( 1 982 )  and 
Brittain  (personal  communication).  Cheatum  (193^)  and 
Horst  and  Costa  (1975)  suggested  that  similar  seasonal 
movements  occurred  in  pulmonate  gastropods. 

Seasonal  changes  in  benthic  distribution  of  Ephemeroptera 
in  lotic  habitats  have  been  reported  by  Ide  (1935) »  Harker 
(1953a),  Macan  (1957),  Ulfstrand  (1968),  Edmunds  et  al . 

(1976)  and  Hall  et  al .  (1980).  In  addition,  Neave  (1930), 
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Haydon  and  Clifford  (197*0,  Keller  (1976),  and  01s  son  and 
Sbderstrom  (1978)  have  observed  substantial  seasonal  upstream 
movements  of  mayfly  larvae  in  running  waters ;  on  the  order  of 
hundreds  of  meters  over  one  or  two  weeks. 

There  is  less  evidence  of  short-term  lateral  movement 
in  lotic  systems.  Elliott  (1971b)  found  that  daily  move¬ 
ments  of  animals  towards  the  centre  of  a  stream  were  more 
pronounced  than  movements  towards  margins.  In  contrast, 

Bird  and  Hynes  (I98I)  reported  that  benthic  movements  of 
lotic  invertebrates  in  an  Ontario  stream  were  random. 

Seasonal  data  on  benthic  distribution  of  E.  inermis 
larvae  in  the  Pembina  River  (Chapter  1)  suggest  that  benthic 
movements  may  indeed  carry  large  numbers  of  individuals 
towards  the  river  centre  at  certain  times.  Maximum  abundance 
of  larvae  occurs  near  shore  during  late  summer,  but  by  the 
time  ice  forms  in  autumn,  densities  have  become  greatest 
in  regions  remote  from  shore.  Because  drift  of  this 
population  is  very  reduced  in  autumn,  this  change  in  micro¬ 
distribution  can  only  occur  by  benthic  movement. 

There  is  no  simple  way  to  estimate  what  degree  of 
lateral  movement  would  be  necessary  to  compensate  for  the 
effect  of  drift  on  distribution  indicated  by  the  model, 
because  there  is  no  unique  set  of  probabilities  of  move¬ 
ment  between  transects  that  will  exactly  counteract  the 
changes  in  benthic  position  of  animals  dictated  by  the 
model's  transition  matrices.  However,  I  used  an  iterative 
method  to  provide  a  general  idea  of  the  minimal  degree  of 
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movement  towards  the  stream  centre  necessary  to  balance 
the  loss  of  insects  from  the  eight  central  transects  of 
the  model  river  after  the  first  day's  drift.  My  assumptions 
were  that  animals  either  remained  in  their  present  transect 
or  moved  towards  the  stream  centre.  The  mean  distance 
moved  was  assumed  greatest  for  animals  in  transects 
adjacent  to  river  margins  (where  current  velocity  was 
minimal)  and  to  decrease  linearly  to  0.1  m  in  transects  12 
and  21  (those  transects  adjacent  to  the  eight  central 
transects) .  No  net  lateral  movement  should  occur  at  the  river 
centre . 

Probabilities  of  movement  of  animals  from  any  particular 
transect  to  other  more  central  transects  were  determined 
according  to  a  Poisson  series.  The  first  iterative  step 
involved  selecting  an  arbitrary  intercept  value  for  the 
relationship  between  initial  distance  from  shore  and  mean 
distance  moved.  The  slope  necessary  to  produce  mean  move¬ 
ment  of  0.1  m  adjacent  to  the  centre  was  then  determined, 
and  the  mean  distance  moved  by  animals  within  all  other 
transects  interpolated  from  this  linear  function. 

Probability  of  lateral  benthic  movement  by  any  animal  from 
a  given  transect  to  each  other  transect  was  calculated 
according  to  the  Poisson  distribution  corresponding  to  the 
interpolated  mean  distance  moved  from  that  transect. 
Probabilities  were  then  assembled  as  a  transition  matrix, 

C' .  The  inner  product  of  this  matrix  and  the  vector 
describing  benthic  distribution  of  animals  after  1  day's 
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drift  (B  )  was  then  calculated, 

B0,  =  B1C '  (27) . 

The  elements  of  B0,  corresponding  to  densities  in  the 
eight  central  transects  of  the  model  river  were  then  summed. 

The  function  describing  benthic  movement  towards  the  river 
centre  was  judged  adequate  to  compensate  for  lateral 
transport  by  drift  if  the  sum  determined  above  was  greater 
than  or  equal  to  the  corresponding  sum  of  elements  of  B°. 

If  the  function  was  inadequate,  a  new  (greater)  intercept 
was  selected,  a  corresponding  slope  determined,  and  a  new 
transition  matrix  calculated  and  evaluated. 

The  maximum  mean  distance  that  would  have  to  be  moved 
by  dead  animals  at  the  river  margin  to  compensate  for 
lateral  transport  by  drift  was  estimated  at  between  0.3  and 
1.0  m  per  day  (Table  16) .  Live  B.  tricaudatus  larvae  would 
have  to  move  a  maximum  mean  distance  of  2-3  m  away  from 
the  river  margin,  depending  on  the  season.  Ephemerella 
inermis  larvae  at  the  riverbanks  would  have  to  travel 
6.6  meters.  Movements  of  this  magnitude  are  probably 
well  within  the  capabilities  of  these  larvae,  provided 
they  are  able  to  orient  in  a  direction  that  will  carry 
them  towards  the  river  centre. 

Elliott  (1971b)  estimated  that  Baetis  rhodani  (Pictet)  and 
Ephemerella  ignita  (Poda)  larvae  moved  modal  distances  of 
6.25  and  4.75  m  upstream  per  day  in  a  trough  placed  within 
a  stream.  Although  Elliott  (1971b)  also  observed  that 
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Table  16.  Functions  describing  estimated  minimum  mean 
lateral  distance  that  must  be  travelled  to 
compensate  for  lateral  transport  by  drift 
in  model  river. 


Species 

Condition 

Intercept 

Slope 

E. 

inermis  (May) 

Live 

6.7 

-0.57 

Dead 

1.0 

-0.08 

B. 

tricaudatus  (May) 

Live 

2.9 

-0.24 

Dead 

0.7 

-0.05 

B. 

tricaudatus  (Oct.) 

Live 

2.1 

-0.17 

Dead 

0.3 

-0.02 

, 


lateral  movements  were  not  as  pronounced  as  were  upstream 
movements ,  he  found  more  animals  crawled  towards  the  stream 
centre  than  towards  shore,  even  though  the  shoreward-facing 
nets  were  situated  very  close  to  the  river  margin.  It  is 
certainly  possible  that  orthokinetic  responses  to  low 
current  velocity  could  result  in  accumulation  of  crawling 
rheophilic  animals  in  faster  water.  Many  studies  have 
indicated  lotic  invertebrates  becoming  substantially  more 
active  when  subjected  to  reduced  current  velocities 
(Edington  19^5>  Elliott  1968a, b).  Animals  in  fast  water 
tend  to  become  oriented  in  an  upstream-facing  direction 
(Neave  1930)  so  that  little  lateral  change  in  position 
should  occur  in  these  regions. 

If,  in  the  absence  of  contrary  evidence,  preferred 
microhabitats  are  defined  as  those  regions  in  which  animals 
are  most  abundant  during  stable  flow  conditions,  then  drift 
does  not  by  itself  result  in  the  maintenance  of  populations 
within  preferred  habitats.  There  is  little  change  in  benthi 
distribution  of  B.  tricaudatus  larvae  during  periods  of 
normal  summer  discharge.  Ephemerella  inermis  larvae  become 
proportionately  more  numerous  near  shore  in  midsummer 
(Chapter  1),  possibly  because  benthic  movements  in  this 
population  do  not  completely  compensate  for  drift.  Thus, 
the  preferred  microhabitat  of  these  populations  does  not 
appear  to  change  dramatically  during  the  summer. 

Nevertheless,  departure  rates  of  larvae,  even  in 
regions  of  maximum  abundance,  are  remarkably  high. 
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Laboratory  experiments  suggest  this  is  unlikely  the  result 
of  density-dependent  interactions.  Departure  from  the  fast- 
water  habitat  may  be  largely  the  result  of  predator  activity 
(Corkum  and  Pointing  1979,  Corkum  and  Clifford  1980, 

Peckarsky  1980b,  Walton  1980b)  or  local  depletion  of 
certain  food  resources  such  as  diatoms  (Hildebrand  1974, 

Keller  1976,  Bohle  1978),  or  both.  Drifting  animals  are  probably 
incapable  of  continuously  monitoring  the  substrate  over 
which  they  pass,  since  the  moment  of. return  to  the  substrate 
is  largely  random  and  dependent  on  chance  contact  with  the 
substrate  (Chapter  3)  •  Thus,  animals  that  settle  in  un¬ 
favourable  regions  must  search  over  the  substrate  to  find 
other  favourable  areas.  Relatively  unidirectional  movements 
(i.e.  towards  the  river  centre)  may  be  a  more  efficient 
means  of  finding  these  microhabitats  than  random  searching 
(Pyke  1978)  ,  especially  if  they  tend  to  run  parallel  to  the 
river  margins. 

The  proposed  system  of  drift  and  lateral  counter¬ 
movement  would  be  particularly  efficient  in  minimizing 
downstream  displacement  of  animals  during  extended  periods 
of  flooding.  During  normal  discharge  conditions,  mean 
downstream  transport  of  all  live  E.  inermis  and  B.  tricaudatus 
larvae  that  drift  are  approximately  1 5  and  8  m, .respectively. 

Net  downstream  displacement  of  B.  tricaudatus  populations 
is  substantially  less,  although  that  of  E.  inermis  is  not. 
Upstream  benthic  movements  may  partially  (Bishop  and  Hynes 
1969b,  Hultin  et  al.  1969,  Elliott  1971b,  Keller  1976, 
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Bird  and  Hynes  1981)  or  completely  (Lehmann  1967,  Marchant 
and  Hynes  I98I)  compensate  for  this.  During  periods  of 
flooding,  however,  current  velocity  at  all  points  in  the 
river  increases  proportionately  to  the  increase  in  discharge. 
Both  the  change  (Butz  1973,  Ciborowski  et  al.  1977)  and  the 
absolute  increase  in  current  velocity  (Madsen  1968, 

Anderson  and  Lehmkuhl  1968)  result  in  increased  ( ’’catastrophic" ) 
drift  and  inhibit  upstream  movement  (Elliott  1971b) . 

Besides  being  more  prone  to  dislodgement ,  animals  will 
require  greater  distances  to  settle.  The  increased 
incidence  of  drift  will  result  in  a  greater  net  transport 
of  individuals  towards  the  river  margins  as  well  as  increased 
downstream  travel.  Thus,  subsequent  erosion  of  animals 
from  their  points  of  return  to  the  substrate  will  produce 
relatively  smaller  downstream  displacement.  As  long  as 
the  erosional  influence  of  the  flood  inhibits  countermove¬ 
ments,  each  dislodgement  will  result  in  further  transport 
closer  to  shore  together  with  a  reduction  in  the  downstream 
distance  travelled.  The  net  effect  on  populations  would  be 
to  produce  greater  abundance  of  individuals  near  shore  and 
reduced  numbers  at  the  stream  centre.  Such  changes  in 
distribution  occur  for  both  populations  during  floods  in 
the  Pembina  River.  Lehmkuhl  and  Anderson  (1972)  observed 
similar  effects  in  an  Oregon  stream  during  flood  events. 

Subsequent  declines  in  discharge  would  again  permit 
access  to  the  river  centre,  and  benthic  movements  would 
result  in  rapid  return  of  individuals  to  ’’preferred"  regions. 
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Drift  and  benthic  movement  are  both  integral  parts  of 
such  a  redistribution  system.  For  this  reason,  one  would 
expect  only  those  organisms  capable  of  free-ranging  benthic 
movement  to  exhibit  a  high  propensity  to  drift.  Indeed, 
relatively  sedentary  taxa  such  as  Chironomidae  and  cyclorrhaphan 
Diptera  do  not,  as  a  rule,  consistently  drift  in  as  great 
proportions  as  do  more  active  animals  such  as  Ephemeroptera 
and  amphipods  (Waters  1972,  Adamus  and  Gaufin  1976). 

The  drift  model  implies  another  ecologically  important 
aspect  of  drift:  the  very  high  turnover  of  individuals  that 
must  occur  throughout  the  river.  If  50  percent  of  the 
animals  depart  from  a  region  each  day,  the  average  residence 
time  of  an  individual  is  2  days.  Benthic  movements  of 
animals  may  produce  either  a  return  of  some  animals  or 
a  loss  of  others.  Townsend  and  Hildrew  (1976)  reported 
a  daily  turnover  rate  of  3*6  percent  for  the  entire  fauna 
of  a  small  stream,  but  as  high  as  43  percent  for  active 
taxa  prone  to  drift.  Sheldon  (1977)  calculated  daily 
turnover  of  Baetis  (mostly  B.  bicaudatus  Dodds)  larvae 
at  approximately  20  percent.  This  compares  favourably 
with  an  estimate  made  by  Elliott  (1971a)  of  the  daily 
proportion  of  Baetis  rhodani  that  drifted  within  a  riffle 
in  a  small  English  stream.  Complete  recolonization  of 
denuded  substrates  by  drift-active  species  in  two  weeks 
or  less  has  been  reported  by  a  number  of  authors  (Waters 
1964,  Ulf strand  et  al .  1974,  Khalaf  and  Tachet  1977>  Gore 
1979,  Wise  and  Molles  1979»  Shaw  and  Minshall  1980) ,  and 


' 

- 

1 

■  S' 


131 


are  indicative  of  daily  turnover  rates  of  at  least  20  percent 
(see  Sheldon  1977)*  Clearly,  such  high  turnover  rates 
indicate  that  individuals  are  capable  of  taking  rapid 
advantage  of  any  favourable  microhabitat  that  becomes 
available . 

Although  the  model's  predictions  of  the  influence  of 
drift  on  distribution  of  animals  pertain  primarily  to  a 
large  river,  the  dynamics  of  transport  and  countermovement 
should  also  be  applicable  to  smaller  streams.  However,  a 
major  difficulty  in  direct  application  of  the.  model  to 
smaller  systems  is  that  mean  current  velocity  is  spacially 
more  variable  in  smaller  streams.  Pools  and  riffles  reflect 
changes  in  lateral  current  velocity  profile  over  relatively 
short  longitudinal  distances;  meanders  in  the  stream  result 
in  changes  in  the  relative  position  of  the  thalweg;  and 
minor  obstructions  such  as  boulders  and  logs  can  profoundly 
interrupt  the  pattern  of  flow.  The  result  is  that  micro¬ 
habitats  available  in  small  streams  will  tend  to  take  on 
a  mosaic  pattern  quite  different  from  the  relatively  homo¬ 
geneous  arrangement  on  which  the  model  was  based.  Any 
animals  transported  directly  downstream  in  the  water  column 
of  a  small  stream  would  be  exposed  to  a  number  of  different 
flow  conditions.  The  instantaneous  settling  rate  would 
still  be  a  function  of  current  velocity,  but  could  vary 
considerably  during  the  course  of  the  drift  episode. 

Estimates  of  departure  from  the  substrate  would  become 
equally  complicated  as  would  lateral  dispersal  determinations. 
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Determination  of  departure  and  settling  probabilities 
would  require  detailed  mapping  of  the  entire  study  area 
and  estimation  of  changes  in  distribution  for  each  point 
in  the  reach  rather  than  for  parallel  transects.  Furthermore, 
an  individual  with  benthic  movement  capabilities  similar 
to  those  postulated  for  animals  in  the  Pembina  River 
would  encounter  a  much  greater  variety  of  microhabitats 
during  one  crawling  episode  than  would  its  counterpart  in 
a  large  river.  Consequently,  one  might  expect  drift  to  be 
of  relatively  less  importance  to  populations  in  small  streams 
than  in  larger  rivers.  Populations  from  large  and  from 
small  lotic  systems  could  possibly  develop  different  behavioural 
tendencies  to  drift  or  to  undertake  benthic  movements. 

A  study  to  test  for  such  differences  would  have  to  be 
carefully  controlled  to  screen  out  confounding  influences 
such  as  variations  in  age  structure  of  the  populations, 
as  well  as  differences  in  the  physical  and  biotic  structure 
of  the  two  systems. 

Drift  of  E.  inermis  and  B.  tricaudatus  larvae  in  the 
Pembina  River  is  of  great  enough  magnitude  to  produce 
profound  changes  in  benthic  microdistribution  of  populations 
over  short  time  intervals.  Short-term  stability  in  distrib¬ 
ution  appears  to  be  maintained  in  spite  of,  rather  than 
because  of,  the  high  propensity  of  animals  to  drift. 

Because  drift  provides  access  largely  to  slower  portions  of 
the  river,  benthic  movement  is  a  necessary  co-requisite  to 
provide  individuals  access  to  all  regions  of  the  lotic  system. 
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Internally  or  externally  mediated  changes  in  the  tendency 
of  individuals  either  to  drift  or  to  engage  in  benthic 
lateral  movements  should  result  in  an  alteration  of  regions 
of  population  concentration  within  the  river.  An  individual 
that  can  rapidly  change  position  in  the  river  minimizes 
downstream  displacement  during  catastrophic  events  while 
at  the  same  time  maximizing  probabilities  of  encountering 
randomly  located  resource-rich  areas.  This  is  clearly 
a  strategy  suited  to  opportunistic  species. 


THESIS  DISCUSSION 


This  study  was  undertaken  to  evaluate  the  role  of 
drift  in  controlling  larval  microdistribution  of  two 
mayfly  populations  in  an  Alberta  river.  I  feel  that 
there  is  strong  evidence  that  drift  is  a  major  determinant 
of  regions  of  greatest  abundance. 

Baetis  tricaudatus  and  Ephemerella  inermis  larvae 
both  exhibit  substantial  levels  of  drift  in  the  Pembina 
River,  yet  distributions  of  these  animals  vary  little 
during  periods  of  stable  flow.  Laboratory  experiments 
suggest  that  the  tendency  to  drift  is  more  strongly  con¬ 
trolled  by  mean  current  velocity  than  by  concentrations 
of  detritus  in  the  substrate  or  benthic  density  of  con- 
specifics.  Live  larvae  are  more  prone  to  appear  in  the 
water  column  than  are  dead  animals ,  suggesting  that 
behaviour  is  important  in  initiating  drift. 

Live  drifting  larvae  rapidly  return  to  the  substrate 
but  tend  to  become  more  strongly  dispersed  laterally  than 
do  dead  larvae.  Because  settling  rates  of  both  live  and 
dead  animals  are  inversely  related  to  current  velocity, 
lateral  transport  during  drift  results  in  animals  being 
most  likely  to  return  to  the  substrate  in  slow  water  near 
shore . 

Although  live  animals  are  more  prone  to  depart  from 
regions  of  slow  current  than  from  faster  water,  this  is 
insufficient  to  counteract  the  effect  of  net  transport 
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towards  shore  during  drift.  This  finding  can  only  he 
reconciled  with  field  observations  on  benthic  distribution 
by  concluding  that  either  the  total  proportion  of  larvae 
that  drift  in  the  Pembina  River  is  very  small,  or  that 
lateral  transport  during  drift  is  compensated  for  by 
benthic  movements  of  the  animals  towards  river  centre. 

I  feel  the  latter  explanation  is  more  plausible  because: 

1)  the  amount  of  drift  normally  occurring  in  the  Pembina 
River  broadly  corresponds  with  or  exceeds  the  levels 

I  would  predict  as  a  result  of  laboratory  and  field 
experiments ; 

2)  flooding,  which  would  be  expected  to  induce  increased 
levels  of  drift  and  to  impede  animals'  benthic  move¬ 
ments,  produces  a  decrease  in  benthic  densities  at 
the  river  centre  and  an  increase  at  the  margins; 

3)  the  only  major  shift  in  population  microdistribution 
observed  in  the  field  during  periods  of  stable  flow 
consisted  of  an  increase  in  abundance  at  stream 
centre  concurrently  with  a  decrease  at  the  margins, 
at  a  time  when  drift  of  animals  within  all  regions 
was  negligible. 

Drifting  animals  with  poor  swimming  abilities  (i.e. 
capabilities  equivalent  to  those  of  dead  organisms  (see 
Elliott  1971a) )  are  transported  relatively  long  distances 
downstream,  but  shoreward  lateral  transport  is  minimized. 

In  contrast,  animals  that  are  capable  of  rapid  return  to 
the  substrate  (e.g.  B.  tricaudatus  and  E.  inermis  larvae) 
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minimize  downstream  displacement  but  are  more  susceptible 
to  transport  towards  stream  margins.  Thus,  rheophilic 
animals  that  drift,  whether  rapid  settlers  or  not,  can 
be  transported  to  regions  of  unfavourable  current  velocity. 

This  is  clearly  of  greater  importance  to  animals  with  poor 
locomotory  abilities,  which  could  not  leave  such  regions 
except  by  further  drift.  Consequently,  one  might  expect 
drift  to  be  less  prevalent  among  relatively  sedentary 
organisms  than  among  individuals  with  well-developed 
crawling  abilities. 

Drift  and  Population  Regulation 

The  prevailing  view  of  the  ecological  significance 
of  drift  is  that  it  serves  as  a  means  of  population  regulation. 
Muller  (195^.  1973,  197*0  and  Waters  ( I961 ,  1965,  1972) 
both  suggest  that  appearance  of  high  numbers  of  organisms 
in  the  drift  results  from  competition  for  food  and/or  space 
in  the  substrate.  Muller  proposes  that  an  important 
auxilliary  function  of  drift  is  dispersal  of  individuals 
to  downstream  reaches  (distributional  drift),  which  requires 
compensatory  upstream  movement  of  adults  prior  to  ovi- 
position.  In  contrast,  Waters  feels  that  drifting  organisms 
simply  represent  excess  production  and  that  most  individuals 
that  drift  could  be  lost  from  the  system.  Both  authors 
view  drift  as  primarily  a  density-dependent  phenomenon, 
but  have  not  elaborated  on  the  nature  of  competition  that 
might  induce  individuals  to  depart  from  the  substrate. 
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There  is  evidence  that  both  interference  and  exploitive 
competition  occur  in  stream  systems  and  can  result  in  drift. 

Several  recent  studies  have  described  aggressive 
interactions  among  space-dependent  tube-dwelling  or 
filter  feeding  organisms.  Hart  (1982)  and  McAuliffe 
(1982)  have  found  that  Hydroptilidae  (Trichoptera)  larvae 
actively  exclude  other  organisms  from  the  immediate 
vicinity  of  their  cases.  Wiley  and  Kohler  (I98I)  and 
Harding  and  Colbo  (I98I)  reported  similar  interactions 
between  black  fly  larvae.  Glass  and  Bovbjerg  (1969) 
and  Hildrew  and  Townsend  (I98O)  presented  evidence  of 
aggressive  behaviour  between  net-spinning  caddisfly  larvae. 
These  interactions  occasionally  result  in  departure  into 
the  water  column  (Wiley  and  Kohler  I98I)  either  as  the 
result  of  active  swimming  movements  or  because  increased 
activity  results  in  erosion  of  individuals  (Hildrew  and 
Townsend  1980)  . 

Steine  (1972)  and  Gyselman  (I98O)  noted  that  mean 
size  of  drifting  Ephemeroptera  was  smaller  than  that  of 
the  benthic  population  and  suggested  this  resulted  from 
competitive  interactions.  However,  the  few  observations 
that  have  been  made  on  behavioural  interactions  between 
mayflies  indicate  that  intraspecific  encounters  produce 
short  distance  crawling  movements  but  not  drift  (Corkum 
1976,  1978a,  Wiley  and  Kohler  1981).  Thus,  interference 
competition  may  or  may  not  be  a  determinant  of  the  magnitude 
of  drift  of  mayfly  populations. 
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My  experiments  on  the  influence  of  benthic  density  on 
departure  rates  of  E.  inermis  and  B.  tricaudatus  larvae 
suggest  that  interference  competition  is  relatively  unimportant 
in  these  populations,  at  least  in  relation  to  the  influence 
of  current  velocity  on  drift.  These  findings  are  consistent 
with  those  of  other  controlled  studies  on  the  relationship 
between  drift  and  benthic  density  (Hildebrand  1974,  Keller 
1976,  Bohle  1978,  Corkum  1978a). 

In  general,  holometabolous  insects,  which  in  many 
cases  have  poor  benthic  locomotory  abilities  and  rely 
on  relatively  restricted  feeding  sites,  may  be  more  prone 
to  drift  as  a  consequence  of  intraspecific  interactions 
than  hemimetabolous  insects  or  crustaceans,  which  sire 
capable  of  free-ranging  foraging  behaviour. 

Drifting  holometabolous  insects  have  poor  settling 
ability  and  generally  behave  similarly  to  dead  animals 
(Elliott  1971a) .  Recurrent  drift  of  these  organisms  will 
result  in  significant  downstream  displacement.  If  settling 
occurs  in  an  unfavourable  region,  the  only  means  of  removal 
is  by  further  drift,  which,  in  the  long  run,  is  counter¬ 
productive.  In  this  sense,  drift  of  these  organisms  is 
consistent  with  Waters'  hypothesis. 

Although  drift  of  highly  motile  animals  may  not  be 
a  primary  response  to  intraspecific  interactions, 
exploitive  competition  could  produce  increased  levels  of 
departure  from  the  substrate.  Drift  often  increases  when 
high  quality  food  (e.g.  periphyton)  is  in  limited  supply 
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(Hildebrand  1972*,  Keller  1976,  Bohle  1978).  Hart  (I98I) 
suggested  that  grazing  pressure  of  the  caddisfly 
Dicosmoecus  gilvipes  Hagen  is  severe  enough  to  make  food 
a  limiting  resource  for  these  animals.  In  such  cases,  drift 
could  regulate  microhabitat  densities  even  though  it  is  not 
a  proximal  response  to  the  presence  of  other  individuals. 

The  frequency  with  which  this  occurs  in  lotic  systems  is 
unknown.  Bishop  and  Hynes  (1969a)  suggested  that  physical 
perturbations  of  the  river  and  predation  may  maintain 
densities  of  animals  below  levels  at  which  food  may  become 
limiting. 

Drift  of  E.  inermi s  and  B.  tricaudatus  larvae  was 
unrelated  to  detrital  concentration  in  my  laboratory  studies, 
even  though  field  population  densities  were  correlated  with 
detritus  levels.  Although  I  did  not  test  for  a  relationship 
between  periphyton  abundance  and  drift,  any  such  relation¬ 
ship  should  not  grossly  affect  the  results  of  the  drift 
model.  Nocturnal  drift  of  both  populations  was  related 
to  mean  current  velocity;  maximal  departure  occurred  at 
lowest  velocities.  Because  periphyton  abundance  in  rivers 
tends  to  be  greater  at  higher  current  velocities  (see 
Horner  and  Welch  I98I),  any  drift  response  to  periphyton 
should  parallel  responses  to  current. 

Drift  and  Foraging 

Both  the  distributional  drift  and  excess  production 
hypotheses  consider  drift  as  a  population  phenomenon  and 
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associate  drift  with  permanent  removal  of  larvae  from 
the  immediate  microhabitat  (Gyselman  I98O) .  This  may 
be  valid  for  sedentary  organisms,  but  my  data  suggest 
that  drifting  larvae  are  not  transported  long  distances, 
and  the  animals  may  be  able  to  maintain  their  position 
in  the  substrate  through  benthic  movements.  If  motile 
organisms  are  capable  of  such  compensatory  benthic  move¬ 
ment,  it  is  unclear  why  they  should  drift.  Drift  is  an 
energy  efficient  means  of  transport  as  compared  to  crawl¬ 
ing;  but  if  an  animal  drifts,  it  must  subsequently  expend 
an  equivalent  amount  of  energy  to  return  to  its  original 
or  to  a  comparable  position  in  the  river.  A  rheophilic 
animal  would  benefit  by  exhibiting  an  orthokinetic  response 
to  current  velocity  whether  or  not  it  regularly  engaged 
in  drift. 

Benthic  movement  is  also  a  more  efficient  means  of 
foraging  over  the  substrate  than  is  drift,  since  a  crawl¬ 
ing  animal  can  continuously  sample  substrate  conditions, 
whereas  a  drifting  animal  can  probably  only  monitor  the 
point  at  which  it  settles.  Clearly,  drift  can  only  be 
advantageous  to  those  animals  capable  of  crawling  when 
benthic  movement  either  becomes  energetically  too  costly 
or  involves  a  greater  risk  of  predation  than  does  drift. 

There  is  strong  evidence  that  crawling  animals  will 
drift  when  in  the  vicinity  of  invertebrate  predators  (Corkum 
and  Pointing  1979 »  Corkum  and  Clifford  I98O,  Peckarsky 
1980b,  Walton  1980b),  yet  E.  inermis  and  B.  tricaudatus 


larvae  frequently  drift  even  when  predators  are  not  present 
One  means  of  examining  the  advantage  of  drift  is  to 
initially  consider  an  active  benthic  animal  that  would 
not  normally  engage  in  drift  behaviour.  Such  an  organism 
should  remain  within  any  microhabitat  that  offers  a 
current  velocity  within  its  zone  of  .preference  and  an 
adequate  food  supply,  and  should  drift  only  if  encountered 
by  a  predator.  A  change  in  current  regime  or  food  avail- 
ibility  will  require  movement  to  a  new  microhabitat. 
Presumably,  this  movement  will  occur  at  night  to  minimize 
exposure  to  visually-feeding  predators  (Allan  1978 » 

Corkum  and  Clifford  1980) .  The  force  of  the  current  will 
tend  to  orient  the  animal  upstream  so  that  most  move¬ 
ments  will  carry  it  against  the  current  (Neave  1930)  and 
possibly  into  faster  flowing  regions,  towards  the  river 
centre.  Movement  should  continue  until  the  animal 
encounters  a  new  microhabitat  with  adequate  resources. 

The  animal  may  be  unsuccessful  in  its  search  for  a  number 
of  reasons: 

1)  favourable  regions  may  be  in  short  supply; 

2)  the  dimensions  of  the  microhabitat  may  exceed  the 
animal’s  moving  abilities; 

3)  there  may  be  a  time  limit  to  the  duration  of  crawl¬ 
ing  movement,  dictated  either  by  the  animal's  energy 
reserves  or  by  the  length  of  the  night; 

4)  the  animal  may  reach  the  upstream  edge  of  the  habitat 
(the  head  of  a  riffle  or  the  margin  of  a  large  stone) 
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A  crawling  animal  that  has  been  unsuccessful  in 
finding  a  suitable  microhabitat  must  either  remain  at  its 
endpoint  (cases  1,  2  and  3)  or  perhaps  retrace  its  path 
(case  4).  Both  possibilities  entail  an  inefficient  use 
of  time  or  energy,  or  both.  Alternatively,  the  organism 
might  drift.  There  may  be  only  a  small  probability  of 
settling  within  a  favourable  area,  but  drifting  would 
require  a  minimum  further  expenditure  of  energy  and  will 
provide  another  opportunity  for  the  animal  to  terminate 
its  search  within  a  suitable  microhabitat.  Whether  the 
final  resting  position  of  the  animal  subsequent  to  drifting 
is  upstream  or  downstream,  and  centreward  or  shoreward,  of 
its  position  prior  to  searching  will  depend  on  the  distance 
previously  moved  over  the  substrate,  the  settling  ability 
of  the  animal  and  on  the  current  velocity  at  the  point  from 
which  drift  occurs. 

Many  authors  have  noted  a  correlation  between  benthic 
activity  and  drift  (Chaston  1968,  Elliott  1968b,  Bishop 
I969,  Waters  1972,  Plosky  and  Brown  I98O,  Hildrew  and 
Townsend  I98O)  and  have  postulated  that  drift  represents 
dislodgement  of  active  animals.  However,  Corkum  (1978b) 
suggested  that  even  though  both  benthic  activity  and  drift 
are  maximal  at  night,  nocturnal  peaks  in  the  two  behaviours 
do  not  generally  correspond;  thus,  drift  and  benthic  activity 
should  possibly  be  regarded  as  independent  phenomena. 

If  drift  does  represent  an  independent  behaviour  that 
follows  unsuccessful  benthic  searching,  one  would  expect 
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nocturnal  drift  maxima  to  lag  behind  peaks  of  benthic 
activity.  The  few  studies  that  compare  drift  and  benthic 
activity  periodicities  suggest  that  drift  peaks  are  equally 
likely  to  precede  as  to  follow  periods  of  greatest  activity 
(Chaston  1968,  Elliott  1968b).  Thus,  although  the  view 
of  drift  as  an  adjunct  to  foraging  behaviour  suggests 
a  mechanism  for  the  initial  evolution  of  active  transport 
in  the  water  column,  its  temporal  relationship  with  benthic 
movement  does  not  appear  to  have  been  maintained. 

Future  Work 

My  interpretation  of  the  influence  of  drift  on 
distribution  of  mayflies  (and  other  motile  taxa)  relies 
heavily  on  the  assumption  that  these  organisms  are 
capable  of  substantial  upstream  and  lateral  benthic  movement. 
Clearly,  this  must  be  adequately  tested  before  the  model 
is  accepted  as  viable. 

Monitoring  of  benthic  movement  of  aquatic  invertebrates 
is  difficult  because  the  organisms  are  small,  cryptic, 
and  predominantly  night-active.  The  flowing-water  medium 
itself  hinders  direct  observation.  Workers  who  have 
attempted  to  study  benthic  movements  by  means  of  traps 
placed  within  rivers  (Bishop  and  Hynes  1969b,  Elliott  1971b, 
Bird  and  Hynes  1981)  have  found  the  movements  to  be  of 
smaller  magnitude  than  drift.  However,  a  shortcoming  of 
these  traps  is  that  animals  crawling  into  them  are  also 
able  to  .crawl  or  drift  out.  More  efficient  trap  design 
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would  alleviate  this  problem  but  would  still  only  produce 
data  on  numbers  of  animals  moving,  not  on  distances  travelled 
by  individuals. 

Other  authors  have  attempted  to  use  mark-recapture 
techniques  to  study  movements  (Harker  1953b,  Brusven  1970, 
Gallepp  and  Hasler  1975,  Hart  and  Resh  1980),  but  such 
methods  are  limited  by  the  potential  effects  of  handling, 
the  large  size  of  populations,  and  an  inability  to  establish 
whether  final  position  has  been  achieved  through  benthic 
movement,  drift  or  some  combination  of  the  two. 

The  use  of  artificial  streams  to  monitor  movements 
of  individual  animals  can  provide  some  estimates  of  move¬ 
ments  (Elliott  1968b,  Corkum  1976),  but  if  animals  do 
indeed  move  to  the  extent  suggested  by  the  drift  model, 
then  most  laboratory  streams  would  be  too  small  to  permit 
individuals  unrestrained  movement. 

Perhaps  the  simplest  and  most  conclusive  test  of  the 
model  and  its  assumptions  would  be  to  impede  drift  or 
benthic  movement  in  the  field  and  to  measure  whether 
predicted  changes  in  distribution  actually  occur.  Benthic 
movement  would  seem  to  be  more  easily  restricted  than  would 
drift.  This  might  be  accomplished  by  partitioning  a  sample 
site  longitudinally  with  slats  embedded  in  the  substrate 
to  restrict  lateral  benthic  movement,  but  still  permit 
lateral  access  to  each  region  by  drift. 

My  model  predicts  that  such  a  partitioning  should 
produce  an  increase  in  densities  near  shore  and  a  decrease 
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at  the  centre,  provided  the  partitions  are  long  enough  to 
prevent  animals  from  crawling  around  the  ends.  This  would 
not  necessarily  establish  whether  or  not  longitudinal 
displacement  of  populations  occurs,  but  it  would  be  a 
measure  of  the  influence  of  lateral  movement  and  transport 
on  maintenance  of  distribution.  By  conducting  the  experiment 
at  sites  with  differing  lateral  slopes,  one  could  also  test 
the  premise  that  increases  in  density  near  shore  should  be 
more  marked  where  depth  and  current  velocity  increase 
rapidly  with  distance  from  the  margin  than  where  there 
is  a  more  gradual  change  in  these  features. 

My  study  suggests  that  E.  inermis  and  B.  tricaudatus 
larvae  in  the  Pembina  River  are  not  static  organisms,  but 
have  the  potential  for  wide-ranging  and  rapid  movement 
throughout  the  lotic  system.  The  voluminous  literature  on 
drift  of  other  taxa  indicates  that  this  is  probably  not 
an  isolated  case.  My  data  also  show  that,  although 
current  velocity  is  important  in  defining  the  microhabitat 
of  these  organisms,  it  is  ultimately  their  behavioural 
responses  to  the  current  which  determine  the  relative 
distribution  of  populations  in  the  river.  This  clearly 
underlines  the  need  for  aquatic  biologists  to  become  more 
aware  of  the  behavioural  capabilities  of  the  organisms  they 
are  studying  to  better  understand  population  dynamics  of 
lotic  taxa. 
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APPENDIX  1 


APL  PROGRAM  TO  DETERMINE  MATRICES  OF  SETTLING 
PROBABILITY  AND  SUMS  OF  PROPORTIONS  DRIFTING 
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[0]  Q+MMATRIX  B\I\ INSERT 

[1]  T  *-Q<-SU  SPENDED+-  (  2  6  26  )p0 

[2]  M+-  (  (  1 1 2  )  -  0 . 5  )  ,  (  8  p  1 2  )  ,  (11  9  7  5  3  1) 

C3]  I+B [1] 

[4]  'GFyFRyll  SETTLING  RATE  INTERCEPT  AND  EXPONENT:' 

[5]  I  NT  1-HU 

[6]  R+INT1 [ 1 ] xM+INTl [ 2 ] 

[7]  'DISPERSION  PARAMETERS: ' 

[8]  T+GENPROB  □ 

C9]  (▼•  •) 

[10]  ( T '  ') 

[11]  (v' RELEASE  POINT  »),(▼!) 

[12]  INSERT+NRNORM 

[13]  Qil\ ]«-26t INSERT 

[14]  SUSPENDED^  I;  ]«-"26  +  INSERT 

[15]  SUSPENDEDII;I1+H  ROU  ND  ( SU  SPENDED  [I;I]+(l-(*(R[I]x"l)  ))+/?[!]) 

[16]  I^J+1 

[17]  +(I<B[ 2])/9 

[18]  Q+EVEN  4  ROUND  Q 


FUNCTION:  Calculates  probability  matrices  K  and  W,  either  for 

entire  model  river  or  for  one  or  more  adjacent  tran¬ 
sects  specified  by  user. 


INPUT : 


OUTPUT: 


Extrinsic:  Right  argument,  a  two  element  vector 

representing  the  first  and  last  transects 
for  which  information  is  desired. 


Intrinsic:  Terminal  will  prompt  for, 

i)  a  two  element  vector  representing 
intercept  and  exponent  values  that 
relate  instantaneous  settling  rate  (r) 
to  distance  from  river  margin  (y). 
These  correspond  to  13 -f  and  131 
(Chapter  4).  ^ 

ii)  a  three  element  vector  representing 
coefficients  that  relate  measure  of 
lateral  diffusion  (s)  to  distance 
from  river  margin  (y).  These  cor¬ 
respond  to  and  13^  (Chapter  4). 

Extrinsic:  A  26x26  element  matrix  corresponding  to 

K.  Rows  for  which  information  was  not 
desired  will  contain  only  zeros. 

Intrinsic:  A  26x26  element  matrix  corresponding  to 

W.  Rows  (transects)  for  which  information 
was  not  desired  will  contain  only  zeros. 
Stored  as  global  variable  SUSPENDED. 


REQUIRED  SUBROUTINES: 


GENPROB,  SD,  PROB,  NRNORM ,  ROUND,  EVEN 


- 


- 


[0]  T+GENPROB  A\I\X 

[1]  7«-(26  2 6  )  p 0 

[2]  I«-l 

[3]  X+-SD  A 

[4]  TliI]+PROB{{I-Q.5) ,XlIl) 

[5]  I«-I  +  l 

[6]  •>(!< 26)/3 


FUNCTION:  Generates  a  26x26  element  matrix,  each  column  of 

which  represents  a  normal  probability  distribution 
with  mean  dictated  by  the  position  of  the  column  in 
the  matrix  and  standard  deviation  externally  specified. 


INPUT: 


OUTPUT: 


Extrinsic : 


Intrinsic : 


Right  argument,  a  3-element  vector, 
entries  of  which  represent  coefficients 
that  relate  lateral  diffusion  measure  (s) 
to  distance  from  margin  of  model  river  (y). 
These  correspond  to  0^,f3^  and  13^  (Chapter  4 

None . 


Extrinsic:  A  26x26  element  matrix.  Each  column  j 

is  a  normal  probability  distribution, 
N(j-0.5»s)  with  entries  (i,j)  the 
probability  of  occurrence  of  a  value 
ranging  between  j-1  and  j. 


Intrinsic:  None. 


SD,  PROB 


REQUIRED  SUBROUTINES: 


' 


[0]  T*SD  A; NET 

[1]  7^26  pO 

[2]  NET*- 1 

[3]  T  lNET]<-Alll+Al2]+AlZ]xMlNETl 

[4]  NET+NET+ 1 

[5]  +(tf£T<26)/3 


FUNCTION:  Calculates  estimates  of  s  for  each  transect  of 

model  river. 


INPUT: 


OUTPUT : 


Extrinsic:  Right  argument,  a  3-element  vector, 

entries  of  which  represent  and  f3^. 

Intrinsic:  None. 

Extrinsic:  A  26-element  vector.  Each  element 

represents  an  estimate  of  s  for  the 
corresponding  transect  of  model  river 
a  distance  of  1  m  downstream  of  a  point 
of  entry  into  water  column. 


Intrinsic:  None. 


REQUIRED  SUBROUTINES: 


None . 


*■  . 


- 
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[0]  L+PROB  B;C;V;A ;F;W 

[1]  L+- 26p0 

[2]  -►(  ( 5  [  2  ]  <  0  )  ,S[2]>0)/3,5 

[3]  L[(B[1]  +  0.5)>1 

[4]  -*0 

[5]  O  0.196854  0.115194  0.000344  0.019527 

[6]  V*-\  2  6-1 

[7]  A+i  (K,K[p7]  +  l)-B[l])xB[2] 

[8]  V+\0 

[9]  F*-A  [  l  +  p7 ] 

[10]  W+1-0 . 5x ( 1++/C* ( |F) *\4 ) *"4 

[11]  V*-V  ,  (  (F<0)  ,F^0)  /  (1-&0  ,W 

[12]  I[l>K[i] 

[13]  -►  (  (p7)=l)/15 

[14]  I[  pK]«*-y[p7]-7[(p7)-l] 

[15]  +((p7)<26)/9 

[16]  I  [  2  6  ]  •♦-I  [26]+(l-V[26]  ) 


FUNCTION:  Generates  a  probability  distribution  bounded  by  values 

0  and  26,  from  a  specified  mean  and  standard  deviation. 

INPUT:  Extrinsic:  Right  argument,  a  two  element  vector 

specifying  mean  and  standard  deviation. 

Intrinsic:  None. 

OUTPUT:  Extrinsic:  A  26-element  vector,  each  element  j ,  of  which 

corresponds  to  probability  of  occurrence  of 
a  value  between  j-1  and  j  according  to 
the  specified  mean  and  standard  deviation. 

Intrinsic:  None. 


REQUIRED  SUBROUTINES:  None. 


. 
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[  0  ]  BOTTOM*NRNORM ;  INTRO  ;  DIST  ;  SPREAD ;  NREM ;  REM ;  BOTTOM ;  ;  SETTLING  ;  I?/?  J 

[1]  DIST* 0 

[ 2 ]  DRIFT*SPREAD*NREM*REM*BOTTOM*26> p  0 

[3]  /?£W[I>1 

[  4  ]  DIST*DIST+1 

[5]  NREM+REM*(*(Rx- 1)) 

[6]  BOTTOM*BOTTOM+  ( REM- NREM) 

[7]  SPREAD+T+ .xNREM 
18]  DRIFT*DRIFT+ SPREAD 
[9]  REM* SPREAD 

[10]  SPREAD* 26 pO 

[11]  -*(  (  +  /REM)>  0.001  )/4 

[12]  BOTTOM *H  ROUND  BOTTOM, DRIFT 

[13]  (j' DISTANCE  TRAVELLED  BY  .001:  • ) , (tDIST) 

[14]  (▼» DISTRIBUTION  OF  SETTLED  ANIMALS:  ’ ) , ( t2 6 + BOTTOM) 

[15]  ( ▼ ’ TOTAL  PROP .  SETTLED:  • ) , (w+ / 2Sf BOTTOM) 

[16]  (  T  ’  TOTAL  PROP .  M777  SETTLED:  ’  ),(t  (  +  /REM)) 

[17]  (t'REL.  PROP .  CAUGHT  IN  A  DRIFT  NET:  ' ),(t(4  DRIF 


FUNCTION:  Calculates  probability  vectors  K.  (probability  that 

an  animal  initially  suspended  in1 transect  i  eventually 
settles  in  transect  j)  and  W.  (sums  of  proportions 
of  animals  originally  in  transect  i  that  will  drift 
past  a  point  on  transect  j  before  settling) .  Also 
determines  downstream  distance  by  which  0.999  of 
animals  have  returned  to  the  substrate. 


INPUT: 


OUTPUT : 


Extrinsic:  None. 

Intrinsic:  Requires  global  variables  R  (vector  of 

instantaneous  settling  rates  of  animals 
in  each  transect;  generated  by  MMATRIX) 
and  T  (matrix  of  normal  probabilities 
generated  by  GENPROB). 


Extrinsic:  A  52-element  vector.  The  first  26 

elements  represent  values  of  K. .  The 
last  26  elements  represent  valies  of 


Intrinsic:  Prints  out  transect  number,  downstream 

distance  by  which  0.999  of  animals  have 
settled,  total  proportion  of  animals  that 
have  and  have  not  settled  at  termination, 
and  vectors  K.  and  W. . 


REQUIRED  SUBROUTINES: 


ROUND 


'  .  5 

' 

v-  ' 
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[0]  A+EVEN  X\R\V\I 

[1]  A+X 

[2]  R+-+/X 

[3]  V+l-R 

[4]  I<- 1 

[5]  i4[J;J]^[J;J]+K[J] 

[6]  I+I+l 

[7]  ^(J<(pK))/5 


FUNCTION: 

INPUT : 


OUTPUT : 


REQUIRED 


Adjusts  diagonal  elements  of  a  square 
that  row  totals  are  1.00. 

Extrinsic:  Right  argument,  matrix  to 

Intrinsic:  None. 

Extrinsic:  Adjusted  matrix. 

Intrinsic:  None. 

SUBROUTINES:  None. 


matrix  such 
be  adjusted. 


[o]  T‘A  ROUND  B 

[1]  T*(10*(-1*A))*|_0. 5+(B“  10*A ) 


FUNCTION:  Rounds  numbers  to  a  given  number  of  places  right  of 

the  decimal. 


INPUT: 


OUTPUT : 


Extrinsic:  Left  argument,  number  of  places  to  right 

of  decimal  desired.  Right  argument,  value(s) 
to  be  rounded. 


Intrinsic : 
Extrinsic : 
Intrinsic : 


None . 

Rounded  value(s). 
None . 


REQUIRED  SUBROUTINES: 


None . 


. 
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APPENDIX  2 


MATRICES  DESCRIBING  SETTLING  PROBABILITIES  AND 
SUMS  OF  PROPORTIONS  OF  DRIFTING  B.  TRICAUDATUS 
AND  E.  INERMIS  LARVAE  IN  MODEL  RIVER 


Each  page  lists  elements  of  two  matrices.  Entries 
in  upper  matrix  describe  probability  that  an  animal  that 
departs  from  substrate  in  any  given  transect  (row)  settles 
in  any  transect  (column).  Entries  in  lower  matrix  describe 
sums  of  proportions  of  all  animals  that  depart  from  a 
given  transect  (row)  that  would  pass  a  point  on  a  given 
transect  (column). 
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irioaudatus  (dead  animals).  Data  for  October  1979. 
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APPENDIX  3 

TEST  RUNS  OF  DRIFT  MODEL  USING  THEORETICAL  DATA 
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Methods 

Behaviour  of  the  model  was  tested  hy  conducting  runs 
in  which  each  component  was  varied  singly  while  holding 
other  components  constant.  The  model  has  three  major 
variables,  each  of  which  are  direct  functions  of  current 
velocity,  thus  indirect  functions  of  distance  from  shore. 
Settling  rate  is  negatively  related  and  lateral  dispersal 
(diffusion)  is  positively  related  to  distance  from  shore. 
Departure  can  vary  positively,  negatively  or  independently 
of  distance  from  shore. 

Individual  values  of  functions  for  each  transect  were 
determined  using  the  distance  from  shore  of  the  midpoint 
of  each  transect.  Benthic  densities  of  animals  used  in  all 
test  runs  approximated  the  benthic  distribution  of  the  two 
mayfly  populations  studied  in  the  Pembina  River.  Animals 
were  most  abundant  at  the  river  centre  (100  individuals 

m  ) .  Abundance  declined  linearly  to  densities  of  four 

_2 

animals  m  in  the  transects  adjacent  to  river  margins.  All 
test  runs  estimated  distributions  of  animals  after  15  days. 

Departure  (p) 

Rate  of  Departure--  Equal  proportions  of  animals  departed 
from  all  transects.  Departure  levels  were  set  at  one  of 
0.2,  0.4,  0.6,  0.8  or  1.0  per  24  hours.  Settling  rate  (r) 
and  lateral  diffusion  (s)  were  constant  at  all  transects 
with  values  of  0.200  and  0.300  respectively. 


■ 
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Departure  as  a  Function  of  Distance  from  Shore--  Departure 
rates  were  made  either  positive  or  negative  linear  functions 
of  distance  from  shore.  Minimum  value  was  assigned  as  0.04. 
Maxima  in  separate  trials  were  0.20,  0.40,  0.60,  0.80  or 
1.00.  Settling  and  diffusion  values  were  constant  as 
indicated  above. 

Settling  Rate 

Settling  rate  normally  varies  inversely  with  current 
velocity,  hence  distance  from  shore  (y) .  Equations  used 
to  calculate  values  of  r  for  the  model  river  were 

r  =  (south  margin,  y<12) 

r  =  f3^(12)fi2  (river  centre,  y>12) 

r  =  13^(2y)fi2  (north  margin,  y<  6). 

Varying  the  Intercept  (B^)--  The  intercept  of  the  above 
functions  is  a  measure  of  animals'  ability  to  settle  at 
all  transects  across  the  river.  When  Bj  is  large,  animals 
settle  over  short  downstream  distances.  As  B^  decreases, 
animals  settle  at  increasingly  greater  distances  downstream 
from  their  point  of  entry  into  the  water  column.  Intercept 
values  (flj)  were  set  at  either  0.50,  0 • 75 »  1.00,  1.25*  1*75 
or  2.00.  The  exponent  (B^)  was  assigned  a  value  of  -1.00 
for  all  trials.  Lateral  diffusion  (s)  was  constant,  with 
a  value  of  Q.300  for  all  transects.  Departure  rates  were 
either  positive  (range  0.04  -  1.00),  negative  (range  0.04  - 
1.00)  or  independent  (0.50)  functions  of  distance  from  shore. 


' 
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Varying  the  Exponent  The  exponent  of  the  above 

functions  is  a  measure  of  the  rate  of  change  of  settling 
rate  (r)  across  the  river.  When  the  exponent  has  a  large 
negative  value,  animals  at  the  river  centre  travel  much 
greater  distances  downstream  before  returning  to  the  sub¬ 
strate  than  do  animals  at  the  margins.  When  the  exponent 
has  a  small  negative  value,  the  differences  in  downstream 
distance  travelled  from  different  portions  of  the  river 
are  minimized.  Exponent  values  ( B^ )  were  set  at  either 
-0.600,  -0.800,  -1.00,  -1.200,  -1.400  or  -1.600.  The 
intercept  (fl^)  was  assigned  a  value  of  1.00  for  all  trials. 
Lateral  diffusion  values  (s)  were  constant  across  the  river, 
with  a  value  of  0 . 300 .  Departure  rates  were  adjusted  as 
indicated  above. 

Lateral  Diffusion  (s) 

Lateral  diffusion  rate  (s)  varies  with  both  the 
distance  that  suspended  animals  have  been  transported  down¬ 
stream  (x)  and  current  velocity,  hence  distance  from  shore  (y) 

1  •  0 .  S  • 

s  =  b3  +  84X  +  G$y 

Because  the  model  determines  values  of  s  at  1-m  downstream 
intervals,  the  above  equation  reduces  to 

s  =  Qj  +  0^(1)°’^  +  B^y 

=  (Bj  +  )  +  B^y 


' 

' 
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The  series  of  equations  used  to  determine  values  of  s  for 
all  transects  in  the  model  river  was 

s  =  (8^  +  8^)  +  O^y  (south  margin,  y<12) 

s  =  (8^  +  B^)  +  8^(12)  (river  centre,  y>12) 

s  =  ( 8^  +  8^)  +  B^(2y)  (north  margin,  y<  6). 

Varying  the  Intercept  (B^+B^)--  The  intercept  of  these 
functions  is  a  measure  of  the  animals'  tendency  to  become 
laterally  dispersed  at  all  transects.  When  the  intercept 
is  large,  animals  become  widely  dispersed.  As  the  inter¬ 
cept  approaches  zero,  all  animals  are  transported  directly 
downstream,  and  there  is  no  lateral  transport.  Intercept 
values  (B^  +  8^)  were  set  at  either  0.100,  0.200,  0.300, 

0.400,  0.800  or  1.000.  The  slope  (B^)  was  set  at  zero. 

Settling  rate  was  constant,  with  a  value  of  0.200  for  all 
transects.  Departure  rates  were  either  positive  (range 
0.04  -  1.00),  negative  (range  0.04  -  1.00)  or  independent 
(O.50)  functions  of  distance  from  shore. 

Varying  the  Slope  (8^)--  The  slope  of  the  lateral 
diffusion  function  is  a  measure  of  the  rate  of  change  of 
lateral  dispersal  across  the  river.  When  the  slope  (B^) 
is  large,  animals  at  the  river  centre  become  more  widely 
dispersed  than  animals  at  the  river  margins.  As  the  slope 
approaches  zero,  dispersal  at  all  points  in  the  river 
becomes  equivalent. 

The  intercept  and  slope  of  the  functions  were  adjusted 


* 
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to  produce  lateral  diffusion  levels  of  zero  at  the  river  margin 
transects  and  to  allow  maximum  lateral  dispersal  values  of 
0.25»  0.50,  O.75,  1.00  or  1.25.  Settling  rate  (r)  was 
constant  at  all  transects  with  a  value  of  0.200.  Departure 
rates  were  adjusted  as  indicated  above. 

Results 

Influence  of  Departure  Rate 

When  both  settling  rates  and  lateral  diffusion  were 
held  constant,  the  distribution  of  animals  after  1 5  days 
was  determined  by  the  sign  of  the  relationship  between 
departure  and  distance  from  shore  (Fig.  29) .  When  de¬ 
parture  was  a  positive  function  of  distance  from  shore, 
densities  at  the  river  centre  became  reduced  and  densities 
increased  evenly  at  all  other  transects.  When  departure 
was  negatively  related  to  distance  from  shore,  densities 
increased  both  at  the  river  centre  and  in  the  transects 
closest  to  shore.  Densities  near  margins  increased  because 
animals  drifting  to  the  river  edges  were  placed  on  the 
substrate  rather  than  being  reflected  back  to  the  river 
centre . 

Increasing  the  maximum  proportion  of  animals  departing 
resulted  in  more  pronounced  changes  in  distribution  after 
15  days.  However,  the  relationship  was  not  linear.  Each 
0.200  increase  in  maximum  departure  rate  produced  a  pro¬ 
portionately  smaller  change  in  distribution. 


V 

. 

. 

>1 


eg 


i 


c 

<D 


Distance  From  Bank,  m 

Figure  29.  Effect  of  varying  daily  proportion  of  animals  de¬ 
parting,  on  benthic  distribution  after  15  days;  settling 
and  diffusion  parameters  constant.  Top,  departure 
increases  with  increasing  mean  current  velocity;  centre, 
departure  independent  of  mean  current  velocity;  bottom, 
departure  decreases  with  increasing  mean  current  velo¬ 
city.  Lines  most  deviant  in  each  figure  are  those  with 
highest  maximal  departure  rates. 
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Influence  of  Settling  Rate 

The  influence  of  varying  both  the  intercept  and 
the  exponent  of  the  settling  rate  functions  was  to 
accelerate  changes  in  distribution,  which  were  determined 
by  relative  departure  rates.  Small  intercept  Values  (Fig. 

30)  and  large  negative  exponent  values  (Fig.  31)  each 
produced  increasingly  marked  changes  in  the  distributions 
after  15  days. 

Influence  of  Diffusion  Rate 

When  the  relationship  between  lateral  diffusion  and 
distance  from  shore  was  constant,  high  diffusion  rates 
resulted  in  more  rapid  changes  in  distribution  than  did 
low  rates.  The  region  of  concentration  of  animals  was 
dependent  on  the  relationship  between  departure  rate  and 
distance  from  shore  (Fig.  32) •  Linear  increments  in  the 
diffusion  values  produced  exponential  changes  in  the 
distribution  of  animals  after  15  days. 

Increasing  the  slope  of  the  functions  relating 
lateral  diffusion  to  distance  from  shore  produced  increasingly 
marked  differences  in  distribution  after  15  days  (Fig.  33) • 
Animals  became  increasingly  abundant  near  river  margins 
as  slopes  became  greater  regardless  of  the  relationship 
between  departure  and  distance  from  shore,  although  when 
departure  was  negatively  related  to  distance  from  shore, 
some  increase  in  density  at  the  river  centre  was  also 


evident . 


. 
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Figure  30.  Effect  of  varying  intercept  of  settling  rate 

function  on  benthic  distribution  of  animals  after  15 
days;  settling  exponent  and  lateral  diffusion  parameters 
constant.  Top,  centre  and  bottom  figures  as  described  in 
Fig.  29.  Lines  most  deviant  in  each  figure  are  those 
with  lowest  settling  rate  function  intercepts. 
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Figure  31*  Effect  of  varying  exponent  of  settling  rate 
function  on  benthic  distribution  of  animals  after  15 
days;  settling  intercept  and  lateral  diffusion  para¬ 
meters  constant.  Top,  centre  and  bottom  figures  as 
described  in  Fig.  29.  Lines  most  deviant  in  each 
figure  are  those  with  largest  negative  settling  function 
exponent  values. 
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Figure  32.  Effect  of  varying  intercept  of  lateral  diffusion 

function  on  benthic  distribution  of  animals  after  15  days; 
lateral  diffusion  slope  and  settling  rate  parameters 
constant.  Top,  centre  and  bottom  figures  as  described 
in  Fig.  29.  Lines  most  deviant  have  largest  lateral 
diffusion  function  intercepts. 
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Figure  33*  Effect  of  varying  slope  of  lateral  diffusion 

function  on  benthic  distribution  after  15  days;  lateral 
diffusion  intercept  and  settling  rate  parameters  constant. 
Top,  centre  and  bottom  figures  as  described  in  Fig.  29* 
Lines  most  deviant  have  largest  lateral  diffusion  function 
slopes . 
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Synopsis 

Departure,  settling  and  lateral  diffusion  all  appear 
to  be  integral  in  determining  the  benthic  distribution  of 
animals  in  the  model  river.  Whether  or  not  animals  become 
concentrated  at  the  river  centre  depends  on  the  sign  of  the 
relationship  between  departure  rate  and  distance  from  shore. 
When  a  greater  proportion  of  animals  depart  from  shoreward 
areas  than  from  the  centre,  populations  will  increase  at 
the  centre.  When  departure  proportions  are  greatest  at 
the  centre,  populations  will  become  most  concentrated  at 
the  margins.  The  process  of  redistribution  is  accelerated 
when  animals  have  consistently  low  settling  rates  or 
consistently  high  lateral  diffusion  values,  or  both. 

Increasing  the  difference  in  settling  rates  between  river 
margin  and  centre  also  accelerates  the  rate  of  redistribution. 
The  degree  of  change  of  the  lateral  diffusion  values  from 
river  margin  to  centre  is  of  critical  importance  to 
eventual  benthic  distribution.  If  values  of  this  parameter 
are  very  small  at  either  river  margin,  it  is  inevitable 
that  the  corresponding  transect  will  eventually  collect 
and  retain  all  drifting  animals  in  the  river.  This  is 
because  a  small  lateral  dispersal  value  can  produce  an 
absorbing  state  in  the  final  transition  matrix.  Thus, 
once  animals  settle  in  corresponding  transects,  the 
probability  of  their  transport  by  drift  from  these  regions 
is  zero.  In  such  cases,  variations  in  departure  or 
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settling  parameters  can  only  accelerate  or  retard  the 
rate  of  accumulation  in  these  areas. 


